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SECTION 1. 0 


INTRODUCTION 

The output power from silicon solar cells employed in space applications 
degrades with time due to electron and proton bombardment. These charged 
particles displace silicon atoms from their normal lattice sites, resulting 
in the creation of recombination centers which reduce the minority-carrier 
lifetime, and hence the-diffusion length. This diffusion-length degradation 
process is commonly characterized by a damage coefficient Kj^, defined 
through the relationship 

I/' - L() - 2 = K l « . (1) 

In this expression Lq and L are pre- and post-irradiation diffusion length, 
respectively, and 0 is the fluence of the irradiating particles (either elec- 
trons or protons) of energy E. Thus, a relatively large Kj^-value is indi- 
cative of vulnerability to radiation damage. 

Recent analyses performed by Brandhorst* indicate that the use of low- 
resistivity (« 1 ohm-cm) bulk silicon for fabricating solar cells should 
result in devices with a considerably larger conversion efficiency than 
that typical of conventional units. Increased open-circuit voltage is the 
primary reason for such expectations. In predicting the radiation response 
of a solar cell fabricated from low- resistivity silicon, the problem arises 
that electron and proton damage coefficients are not known for such material 
and thus experimental determinations of these quantities are required. 

In this report, we present the results of an experimental program aimed 
at determining electron and proton damage coefficients in low-resistivity 
boron-doped silicon. Minority-carrier lifetimes were measured before 
and after irradiation for bulk silicon samples using a steady-state photo- 
conductivity technique. In addition to measurements on 0. 1 ohm-cm 
material, comparison data were obtained for higher resistivity samples, 
thus permitting determination of the dependence of damage coefficient on 
resistivity. Pre- and post-irradiation measurements of two types were 
also made for n-on-p solar cells of various base resistivities, including 
0. 1 ohm-cm: 1) diffusion length measurements using penetrating radiation 

(Co-60 source); 2) determination of cur rent- voltage characteristics during 
illumination of cells with a solar simulator. For low- resistivity bulk 
silicon, the possible dependence of on growth technique (oxygen con- 
centration) and dislocation density was examined. Electron irradiations 
were performed at three energies (0. 5, 1.0, and 2. 5 MeV) to permit 
examination of the energy dependence of K-^. Following irradiation, 
isothermal annealing studies were conducted at room temperature and 
60°C for selected bulk samples. In addition to experimental findings, 
analytical results pertaining to certain aspects of the investigation are 
also presented. 
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SECTION 2. 0 


EXPERIMENTAL TECHNIQUES 

In this section, experimental techniques employed and additional related 
information are presented. Included are descriptions of: a) measurement 

techniques; b) sample preparation procedures; c) radiation facilities and 
irradiation procedures. 

2. 1 Measurement Techniques 

The ASTM standard method for measuring minority-carrier lifetime in 
bulk silicon is the photoconductivity decay technique. ^ However, imple- 
mentation of this technique becomes increasingly difficult as resistivity 
is decreased. Furthermore, data interpretation can be a problem when 
trapping effects are significant. In an attempt to overcome such problems, 
we have devised a steady- state method of lifetime measurement for low- 
resistivity material. As shown in Figure 1, a bulk silicon specimen with 
ohmic contacts attached is situated in an enclosed chamber and is illumi- 
nated by penetrating light from a steady-state source (projection lamp- - 
GE DCA). Non-penetrating wavelengths are eliminated by a silicon filter 
considerably thicker than the test sample, and therefore a condition of 
uniform carrier excitation is obtained. A Corning 7-57 filter was em- 
ployed to reduce heating of the Si filter. A variable-frequency light chopper 
(0-300 Hz) is employed between the specimen and the projection lamp. A 
chopping frequency of 277 Hz was utilized for all lifetime measurements 
in this program. The purpose of chopping the light is to permit the use of 
phase - sensitive detection - a convenient approach to enhancing the signal- 
to-noise ratio. The light chopper also provides the reference signal in 
such a detection scheme. Also shown in Figure 1 is a silicon photodiode 
which provides a signal proportional to light intensity. Background lamps 
(unchopped) were employed to fill traps, and this aspect of the experiment 
is discussed in Section 3. 7. In order to avoid generation of a contact 
photovoltage, sample contacts were shielded from the chopped light, as 
shown in the figure. Approximately the central third of a sample was 
illuminated. 

Figure 2 shows a schematic of the electrical system employed. The sample- 
under-test is connected in a constant-current circuit (R »R g ), and the steady- 
state photoconductivity (SSPC) signal, AV, due to chopped-light illumination 
is monitored. Voltage gain is provided by a low-noise transformer (Triad 
Geoformer G-10) and an operational amplifier (Analog Devices 118A). The 
signal is phase- sensitively detected by means of an Ithaco Model 353. Using 
appropriate circuit considerations, the measured signal was related to the 
signal of interest, AV. 
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Figure 1. Schematic representation of the steady-state photoconductivity (SSPC) apparatus 
employed to measure minority- carrier lifetime in bulk silicon specimens. 
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SSPC apparatus to determine minority- carrier lifetime in bulk silicon specimens. 




For a constant current circuit, excess minority-carrier density An (assum- 
ing p-type material) is given by the relation 


An 


AV P o^h 
V (U + pJ 

o e h 


AV 

eV o p(n e +U. h ) ' 


( 2 ) 


where p Q is equilibrium majority carrier density, p is resistivity, V Q is the 
dc voltage drop across the illuminated portion of the sample, e is the elec- 
tronic charge, and p^ and p g are hole and electron drift mobilities, respec- 
tively. (Equation 2 assumes that AV << V , a condition that held true for 
the present study. ) Measured excess carrier density from an SSPC experi- 
ment is related to minority- carrier lifetime t by 

An = gT, (3) 

where g is carrier generation rate. The generation rate for the SSPC ap- 
paratus was determined using a calibration procedure described in Section 
3.6. We measured An and g was known, thus yielding t. To obtain diffusion" 
length L, the expression 

L = (Dt) 1/2 (4) 

was employed. The procedure used to determine diffusivity D is dis- 
cussed in Section 3. 2. All lifetime measurements using the SSPC apparatus 
were made at room temperature (2 5 ±5°C). 

Measurements of lifetime for several unirradiated specimens were also made 
using transient photoconductivity decay techniques. A flash x-ray source 
was used for carrier excitation; the experimental apparatus employed has 
been described previously. ^ 


Measurements of diffusion length were made for n-on-p Si solar cells using 
the Co-60 gamma-ray source at Northrop Research and Technology Center. 
This measurement method has been described by Gremmelmaier, ^ Rosen- 
zweig, ^ and, more recently, by Reynolds and Meulenberg. ^ Measurements 
of short-circuit current I sc were made and then diffusion lengths were de- 
termined using the expression6 


I 

sc 


q A g Y L 


cosh(w/L) - 
sinh(w/L) 


i 


(5) 


where A is cell area, w is cell width, and g^ is the carrier generation rate 
for the Co-60 source. This rate was determined from a series of total 
dose measurements using CaF thermoluminescent dosimeters. Measured 
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dose rate in rads(Si)/sec was converted to using the conversion factor 
4. 2 x 10 13 pairs /cm 3 rad. The result obtained was g^=6.55x 10^ pairs/ 
cm 3 sec. 

Equation (5) only holds for the case of an infinite surface recombination 
velocity (s) at the back surface of a solar cell. We assumed infinite s 
here. Using a more general expression^ to calculate L for various values 
of s reveals that the difference between results for s = 10^ cm/sec and 
s = oo is very small. ' r The assumption of infinite s at an ohmic contact 
seems quite reasonable. We finally note that employment of Equation 5 
to determine L based on measured I__ requires an iterative procedure, 
and this was accomplished using a programmable desk-top calculator. 

Current- voltage characteristics for solar cells were obtained before and 
after irradiation using the solar simulator at Centralab Semiconductor. 
Measurements were made under AMO conditions at room temperature 
(2 5 to 30°C). 

An alternative measurement technique that was considered for employment 
here is the steady- state surface photovoltage (SPV) method. 7- 10 -phis 
approach to measuring diffusion length was developed quite a few years 
ago, but only recently has it begun to receive attention by various research 
groups. In the present investigation, a comparison was made for several 
samples of diffusion lengths determined using the SSPC method and the 
SPV method. ** v 

2. 2 Sample Preparation Procedures 

Low- resistivity (' v 0. 1 ohm-cm) boron-doped silicon ingots for use in this pro- 
gram were obtained from Dow Corning. Four ingots were obtained, corres- 
ponding to the four possible combinations of high and low oxygen concentration 
and high and low dislocation density. Two ingots were grown by the float- zone 
method and were expected to have a relatively low oxygen content. The other 
two ingots were grown by the Czochralski method (relatively high oxygen con- 
tent). One float-zone crystal and one Czochralski crystal were grown to be 


For the interested reader, we note that the general expression referred 
to here (Equation A9 in Reference 6) appears to have a few typographical 
errors. It should read (using the present notation): 

I = a A p T T Pexpt-.w/LU: £ L |exp(-w/L) ^l[ 
sc ®y |_ s L sinh(w/L) - D cosh(w/L) 



SPV measurements were made at the Jet Propulsion Laboratory through 
the cooperation of R. J. Stirn. 
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dislocation free, and the other two ingots had a high dislocation density. 
(Measurements of dislocation density are described in Section 3. 4). Studies 
were also made on higher resistivity material, also grown by Dow Corning. 
This material was already on hand at Northrop Research and Technology 
Center at the start of this program. 

An important factor in deciding on bulk- silicon test-specimen dimensions 
was the bombarding particle energies to be employed. Table I shows the 
energies used and the corresponding particle ranges in silicon. ^ It is 
desirable to avoid a condition of inhomogeneous damage for the purpose 
of aiding data interpretation. We selected a sample thickness of 0. 3 cm, 
which is comparable to that employed in conventional solar cells. This 
thickness is considerably less than the particle ranges under consideration, 
and thus reasonably uniform damage profiles should result for 1. 0- and 2. 5- 
MeV electrons and for 10-MeV protons. For 0. 5-MeV electrons, although 
the range is three times larger than the sample thickness, inhomogeneities 
are expected to occur because the variation of damage coefficient with particle 
energy is quite strong near 0. 5 MeV. Upon penetrating half the thickness of 
a 0. 3 mm sample, a 0. 5-MeV electron loses ~0. 057 MeV. The 0. 5-MeV- 
to-0.443-MeV damage coefficient ratio for silicon^ is ~ 1.45, and thus the 
amount of damage expected at the center of the sample is ^70% of that at 
the bombarded surface. (A similar calculation for incident 1. 0-MeV elec- 
trons yields ~90% at the sample center. ) 


TABLE I 


Electron 

Energy 

(MeV) 

Electron 

Range 

(mm) 

Proton 

Energy 

(MeV) 

Proton 

Range 

(mm) 

0. 5 

0. 94 

10 

0. 71 

1 . 0 

2. 3 



2. 5 

6. 2 




Electron and proton ranges in silicon for the 
bombardment energies employed in this study. 


Bulk silicon samples measuring approximately 1 cm x 2 cm x 0. 03 cm 
were fabricated at Centralab from several Northrop- supplied silicon 
ingots. Centralab also prepared samples for our use from 6 ohm-cm 
Czochralski material (B-doped) that is conventionally employed in the 
fabrication of solar cells. Additionally, n-on-p solar cells were fabricated 
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at Centralab using the same Northrop- supplied ingots employed in the pre- 
paration of bulk samples, thus permitting a meaningful device-bulk com- 
parison. Conventional production-line n-on-p cells (Czochralski) were 
also obtained from Centralab for the purpose of providing a "standard- 
cell 11 comparison during the radiation testing. 

All bulk- samples were sandblasted at Centralab for the purpose of achieving 
a high surface recombination velocity. (Surface recombination considerations 
are discussed in subsequent sections of this report. ) Ohmic contacts were 
prepared at Northrop using ultrasonic soldering techniques. Contacts (two 
per sample) were narrow (^0. 1 cm) stripes prepared on a major face at 
the ends along the 1-cm dimension. Contacts were checked by making 
forward and reverse resistance measurements at room temperature. Re- 
sistivity was measured for each specimen using conventional four-point- 
probe techniques. 

2. 3 Radiation Facilities and Irradiation Procedures 

Electron irradiations were performed at 0. 5“-, 1. 0-, and 2. 5-MeV using 
the Dynamitron at the Jet Propulsion Laboratory (JPL), Pasadena, Calif. 

All bulk specimens and solar cells were illuminated during irradiation at 
an intensity of 140 mW/cm^ using the solar simulator at the Dynamitron 
facility. Test specimens were bombarded while in vacuum at a pressure 
of ^10”^ mm Hg. Samples were attached to an A1 plate using Apiezon H 
vacuum grease to provide a good thermal contact (see grease discussion 
in Section 3. 5). The temperature of the plate was controlled at 25°C. To 
minimize room-temperature annealing, all irradiated bulk samples and 
solar cells were stored at dry ice temperature following irradiation, with 
storage beginning approximately one-half hour after bombardment. " Spec- 
imens were removed from dry ice for about an additional one-half hour 
to permit a measurement, then were returned to dry ice until time for the 
next irradiation. Current-voltage characteristics for solar cells and life- 
time measurements for bulk samples were obtained after a total of ~ 1 hour 
at room temperature. Diffusion length measurements on selected solar 
cells using a Co-60 source were made after ^2 hours at room temperature. 

At the Dynamitron, a Faraday cup was employed to determine the fluence. 

The following fluences were utilized at all three electron energies: 3 x 10 1^, 
10^4, 3 x 10^, 1015, and 3 x 10 1 5 electrons /cm^. For the first three 
fluences, the flux was 5 x 10^0 electrons /cm^ sec. At 10 1 ^ electrons/ 
cm^, the flux was loH electrons /cm^ sec; and at 3 x 10^ electrons /cm^, 
it was 2 x 10 11 electrons /cm^ sec. The flux was increased with increasing 
fluence in order to obtain reasonable exposure times. (The last exposure 
at a given energy required 2. 8 h. ) 
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The 13-MeV Tandem Van de Graaff at the California Institute of Technology, 
Pasadena, California, was employed for 10-MeV proton irradiations. Bulk 
samples and solar cells were illuminated during irradiation at an intensity 
of 140 mW/cm^ using a Varian VIX-150 Xenon Short- Arc Illuminator. A 
Coherent Radiation Model 201 Broad Band Power Meter was employed to 
calibrate the xenon illuminator and map the light-beam profile. At the 
proton source, because of beam uniformity and beam size considerations 
only one sample could be irradiated at a time. The xenon lamp could be 
employed to illuminate a single sample to an intensity of 140 mW/cm^ to 
within about ±5% over the specimen surface. Samples were irradiated in 
vacuum at a pressure of ~10 - ^ mm Hg. The light source was external to 
the vacuum chamber and illumination was accomplished through a fused 
quartz window (Type 12 5) in the chamber wall. This window has a trans- 
mittance of nearly 100% over the wavelength range from 0. 2 7 to 2. 7 pm. 
Samples were mounted with Apiezon H vacuum grease on a massive, moveable 
block of A1 in the chamber and externally moved into the beam one at a time. 
Irradiation times were ^30 seconds per sample for all fluences. The tem- 
perature rise of the aluminum block was determined to be less than 1°C 
during the irradiation of each lot of about 15 samples. Furthermore, the 
temperature rise of the sample relative to the sample mount was determined 
to be about 1°C when subjected to the light. In the absence of temperature 
control mechanisms, the sample holder was, therefore, brought to an initial 
temperature of 2 3°C prior to evacuation, giving sample temperatures ranging 
from 24 to 25°C during irradiation. Fluences employed were the following: 

~3 x 1010, ~10", ^3 x 10", and 1.2 x 10" p/cm^. 

Fluences were controlled by means of a current integrator which monitored 
the charge impinging on the sample mount through an aperture of known area. 
The sample holder was maintained at positive bias to collect secondary elec- 
trons. Integrator offset current and leakage current due to finite resistance 
from sample to ground were balanced out, and the beam currents employed 
were chosen to be large compared to residual drift in these quantities. A 
thin wire was scanned past the aperture to measure beam current as a function 
of position and verify beam uniformity over the aperture. This method of 
regulating fluence was deemed to be preferable to the use of a Faraday cup 
since the area of beam uniformity was not much larger than the sample size. 

A Faraday cup would have had to be removed from the area of interest, and 
would have had a very small beam aperture, necessitating the measurement 
of much smaller currents. The fluxes for proton irradiations were as follows: 
The fluence of 3 x 10 1 0 p/cm^ was obtained at a flux of about 6 x 10^ p/cm^ 
sec; the fluence of 10^ at 6-9 x 10^; of 3 x 10" at 2 x 10^0; and the fluence 
of 1.2 x 10" p/cm^ at a rate of 3. 6 x 10 ^ p/cm^ sec. The fluxes were 
chosen so as to keep total time at room temperature within acceptable limits 
of about 1 h for each group of samples which were mounted in vacuum simul- 
taneously. Dry ice storage was utilized in the same manner as described 
above for electron irradiations. 
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2. 4 Sample Designation System 

As discussed above, a set of bulk samples and solar cells were fabricated 
from the same silicon ingot to permit a meaningful comparison between 
device and bulk behavior. This same procedure was followed for a variety 
of silicon ingots. The sample designation system employed in this report, 
both in the text and on all figures, provides information concerning growth 
technique, resistivity, dislocation density, and whether a specimen is a 
bulk sample or a solar cell. For float-zone material, the abbreviation FZ 
is employed, and for Czochralski material CZ is used. We employ the 
shorthand designations HD for samples prepared from material with a high 
dislocation density and LD for low-dislocation-density units. . (These des- 
ignations only appear for low-resistivity specimens because no dislocation- 
density determinations were made for higher-resistivity material. ) Nom- 
inal resistivities are given in ohm- cm and also a number is assigned to 
each bulk sample and device. 

A few examples will serve to illustrate the designation system: a) FZLD 

0. 1/4: this notation is for a bulk silicon sample (number 4) fabricated from 
0. 1 ohm-cm (nominal) low-dislocation-density float-zone material; b) 
CZHD 0. 1/SC2: this designation is for a solar cell (SC), device number 2, 
fabricated from 0. 1 ohm-cm high-dislocation-density Czochralski material; 
c) FZ 13/4: a bulk sample (number 4) fabricated from 13 ohm-cm float- 
zone material. Bulk samples designated CZ 6/x were fabricated from 
material typically employed at Centralab for the fabrication of solar cells. 
Solar cells designated CZ 9/SCx are typical production line cells obtained 
from Centralab. 

Data points on all figures are identified in two ways. Some symbols are 
associated with measurements on a particular bulk sample or solar cell, 
and in these cases the sample (or device) number is listed on the figure 
in a legend. Other symbols represent an average of data obtained for two 
or more bulk specimens of the same type, and in these cases the legend 
identifying this type of data point will list the ingot designation only, with 
no specimen number appearing. For example, a data point identified as 
CZLD 0. 1 represents an average value for identically treated bulk samples 
from ingot CZLD 0.1. 
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SECTION 3. 0 


PRELIMINARY EXPERIMENTS AND ANALYSES 

In the course of preparing for irradiations and subsequent determination of 
damage coefficients, several preliminary experiments were performed. Ad- 
ditionally, analyses were made of certain features of the experiment, includ- 
ing surface recombination effects and other phenomena of interest. This 
section discusses the various preliminary considerations that were made. 

3. 1 Analysis of Surface Recombination Effects 

The employment in this investigation of thin bulk silicon samples with sand- 
blasted surfaces necessitates accounting for surface recombination effects. 
(Justification for choosing sandblasted surfaces, along with a description 
of associated experimental results, is given in Section 3. 8. ) In general, 
one measures an effective lifetime t 6 ££ which may or may not be equal to 
the bulk lifetime t depending on the degree of surface recombination occur- 
ring. An analysis of such effects was performed for the sample geometry 
used in this study. The sample can be approximated by a thin slab having 
one finite and two infinite dimensions. The problem is to determine the 
ratio of T e ff to t for a given value of bulk diffusion length L and a given 
surface recombination velocity s. 

For a thin slab of thickness 2T under conditions of uniform carrier excitation, 
the following expression holds: LJ 


T elf /T = 


(sL/DT) sinh (T/L) 

(1/L) sinh(T/L) + (s /D) cosh (T/L) * 


( 6 ) 


Calculations appropriate for 0. 1 ohm-cm p-type s^licon^were made using 
Equation 6 by assuming a diffusivity D of 11. 7 cm /sec and a sample half- 
thickness of 165 (im. Computational results are presented in Figure 3, 
where T e ff/ T is plotted versus bulk diffusion length with surface recombina- 
tion velocity as a parameter. Assuming that the pre-irradiation value of 
L is 100 urn, for the largest values of s measured lifetime will differ from 
the actual value by as much as a factor of about two. On the other hand, 
if s is ^ 100 then the error will be no larger than about 5% before irradiation. 
Another noteworthy feature of the calculations is that the difference between 
results for s = 10^ cm/sec and s = 00 is very small. Thus, as long as the 
condition s >10^ cm/sec holds for specimen samples, differences in s 
from sample to sample with otherwise identical properties will be only of 
second-order importance in terms of measured lifetime ( T e 

' n This value of D is employed only to illustrate the calculational procedure. 
As noted in Section 3 0 2, a slightly smaller electron diffusivity was utilized 
in this investigation for all calculations relating to experimental findings. 
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Figure 3. The ratio of effective to actual lifetime versus diffusion length with surface recombination velocity 
s as a parameter for a thin silicon slab. The solid curves were calculated using the relationship 
shown along with the listed values of diffusivity D and half- thickness T. The dashed curves were 
calculated using the two approximations shown (see text for discussion). 



Since the variation of T e ff/T is approximately independent of surface recom- 
bination velocity for s ^ 10^, analysis is simplified for samples with such sur- 
face properties. Two examples will illustrate this point. First, we note 
that for T > L, sinh (T/L) « cosh (T/L). If it is further assumed that sL 
>> D, Eq. (6) reduces to 


eff 


h 



(7) 


This expression was used to obtain one of the dashed curves in Figure 3, 
assuming T = 1 6 5 Mm, and agreement with the s = 05 curve is seen to be 
excellent for L ^ 100 Mm. The deviation for larger L is due to the deviation 
of sinh (T/L) from cosh (T/L). It is thus seen that for the case of large s 
the analytical dependence of T e ff/ T on diffusion length assumes a very simple 
form. 


From the standpoint of correcting for surface recombination effects, it is 
desirable to express t as 



( 8 ) 


where T g is "surface" lifetime. Such an expression implicitly assumes a 
constant value for T g . It turns out that this approach is strictly true only 
' for the case of T « L. For this situation, Eq. (6) yields t = T/s. For 
our case the inequality T « L will certainly not hold. However, Eq. (8) 
constitutes a very good approximation for the present situation, as seen 
by the following example. 

Consider a sandblasted 0. 1 ohm- cm' p- type sample with a half -thickness T 
of 165 am, a pre-irradiation diffusion length of 100 pm, and a surface re- 
combination velocity of 10^ cm/sec. Assuming a diffusivity D of 11. 7 cm^ / 
sec, we obtain a pre-irradiation value for t of 8. 5 psec. From Figure 3, 
t ff /T = 0. 49 for L = 100 pm, and therefore T e ££ = 4. 2 psec. Equation (8) 
yields T g ~ 8. 2 psec. We next assume that T g remains constant as T de- 
grades due to irradiation. The second dashed curve in Figure 3 shows 
the results of calculating t^/t for various L. values (<100 pm, normalized 
to solid s = 10-4 curve at 100 am) by employing Eq. (8) and assuming T g - 
8. 2 psec. It is seen that the largest deviation of the dashed (approximate) 
curve from the actual (exact) curve is only on the order of 5%. Therefore, 
one could accurately employ Eq. (8) to interpret steady- state photoconductivity 
measurements. The same surface correction that is applied prior to ir- 
radiation can be applied following bombardment. 
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Under the assumption of a constant surface lifetime, T , it is readily shown 
that if effective values of diffusion length are employed to determine damage 
coefficients, then Equation 1 takes the form 


(L 


eff -2 


(L 


eff -2 
0 ' 


K 0 . 
L 


(9) 


* GXX 

If we denote a damage coefficient determined with Equation 9 as , then 
based on the above discussion one would expect K^ e 1 to be in good agree- 
ment with a damage coefficient determined using actual bulk diffusion lengths 
(Kj^). This was found to be the case experimentally. Values for were, 

on the average, ^6% larger than those for (based on data obtained follow- 
ing 1-MeV electron irradiation). 


Although we could have employed the above-described effective-diffusion- 
length approach to determining in the present investigation, the more 
informative approach of obtaining actual pre- and post-irradiation diffusion 
lengths was followed. This required that each measured (effective) lifetime 
be converted to actual bulk diffusion length by properly correcting for the 
effects of surface recombination. For this purpose, Equation 6 was em- 
ployed in a slightly different form: 


T eff D 

Appropriate values of s and D (both discussed below) were employed and that 
value of L required to satisfy Equation 10 was then determined by using an 
iterative procedure on a programmable desk-top calculator. 


1 - 


(sL/DT) sinh (T/L) 


(1/L) sinh (T/L) + (s/D) cosh (T/L) 


(10) 


3. Z Procedure for Determining Diffusivity and Mobility 

Determination of diffusion-length damage coefficients was a primary goal 
of the present study. Because minority-carrier lifetime is the measured 
quantity, utilization of Equation 4 is required to determine diffusion length, 
which requires knowledge of diffusivity D. Additionally, determination of 
lifetime with the SSPC method requires knowledge of drift mobilities. (Of 
course, u, and D are related through the Einstein relation. ) Since Lifetime 
measurements were made as a function of resistivity, and thus D varied, 
it became necessary to assume reasonable values for diffusivity in a con- 
sistent manner (i. e. , consistent with measured resistivities). The follow- 
ing approach was employed. Irvin’s data^ was used to determine acceptor 
concentration N a from measured resistivity. Next, hole drift mobility 
was calculated using the expression 

u = 1/oeN (11) 

h a 


14 



which assumes N a = p Q . Next, Irvin's data was again employed to determine 
the resistivity of n-type Si with the same impurity concentration as determined 
above for p-type material. Equation 11 then yields for the new value of p 
(with donor concentration substituted for N a ). This approach cannot be 
directly employed for 0. 1 ohm- cm material because all impurity atoms are 
not ionized. We determined the number of ionized impurities, based on the . 
value of N determined from Irvin, in the following manner. 

For 0. 1 ohm-cm boron-doped silicon, the equilibrium Fermi level Ej is suf- 
ficiently close to the acceptor level E a so that a significant fraction of the 
acceptor atoms are not ionized at room temperature. The ratio of ionized 
acceptor concentration N a to the total acceptor concentration N a is given 
by 15 


N ’ i 

a_ _ 1 n? , 

N ~ 1 + 3 exp [ (E - Ej/kT] ’ 1 ' 

a a I 

where 3 is the impurity level spin degeneracy (3=4 for acceptor levels^). 
For 0. 1 ohm-cm material, N a is about 4. 8 x 10^ cm’^ (from Irvin^). 
Additionally, E a is located 0. 045 eV above the valence band. By consider- 
ing the variation of Ej with carrier concentration simultaneously with Equa- 
tion 12, the following results were obtained: E£ is located approximately 
0. 092 eV above the valence band and about 60% of the acceptor atoms are 
ionized. 

Once the ionized acceptor density is determined by this method, then Equa- 
tion 11 yields (substitute N a for N a ). Following this approach yields 

M, - 210 cm^/V-sec for 0. 1 ohm-cm material. Determination of Mu for 
n e 

low-resistivity silicon using the above approach presents somewhat of a 

problem. One would like to determine the electron mobility in p-type 
material with N acceptor atoms, N~ of which are ionized. Irvin's plot 
of p vs N for n- and p-type Si has built into it certain assumptions regard- 
ing impurity ionization for low-resistivity material (spin degeneracy, im- 
purity level locations, etc. ) Without knowledge of these precise assump- 
tions, one hesitates to attempt a detailed calculation for physically equi- 
valent n-type material (equivalent in terms of scattering effectiveness) 
which would yield, in conjunction with Irvin's graphs, a value for u< e in 
0. 1 ohm-cm material. Instead, we have employed the approach of aver- 
aging literature values for this quantity. Sze and Irvin^° give a value of 
~435 cm^/V-sec and Wolfstirn's value^ is ~405 cm^/V-sec. The av- 
erage of these two quantities is numerically equal to that given in Neu- 
berger^, ~420 cnW/V-sec, and hence we employ here this value for p e 
in 0. 1 ohm-cm p-type Si. The corresponding electron diffusivity is 
10. 9 cm^/sec at 300°K. 
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It may not be clear yet to the reader why we employed the approach described 
above for obtaining mobility values in silicon of various resistivities. To 
determine An using Equation 2 mobility values must be assumed. Commonly 
employed references, such as Sze and Irvin^, give graphs of mobility vs 
impurity concentration N. Thus, to use such graphs one must determine N. 

A common approach is to measure resistivity and then use Irvin's graph^ 
to determine N, followed by the obtaining of mobility values from Sze and 
Irvin. However, built into Irvin's graph are assumed mobility values which 
are not completely consistent with the mobility graphs of Sze and Irvin. (In 
fact, there is an apparent plotting error in Sze and Irvin's curve in that 
at low impurity concentrations a limiting value of ~600 cm^/V-sec is ap- 
proached, as opposed to the accepted value of ~500 cm^/V-sec. ) Hence, 
for purposes of consistency we decided to employ Irvin's graph alone to 
determine mobility, as described above, with the exception of the u e value 
for 0. 1 ohm-cm material. (This approach removed a considerable amount 
of scatter initially observed in a comparison of SSPC lifetimes and photo- 
conductivity decay (PCD) lifetimes performed during the generation rate 
calibration (described below) for the SSPC apparatus. ) 


As an experimental check on the hole mobility assumed for 0. 1 ohm-cm ma- 
terial, Hall mobility measurements were made at room temperature using 
conventional techniques for one sample from each of the four low- resistivity 
ingots being studied in this program. Results are shown in Table II, along 
with calculated drift mobilities. (A Hall-to-drift mobility ratio 0. 92 was 
assumed. ^®) Agreement among the four ingots is good. It is satisfying to 
note that the average hole drift mobility in Table II (223 cm^/V sec) is in 
good agreement (~6% higher) with the value of 210 cm^/V sec determined 
above. 

TABLE II 


Ingot 

FZLD 0. 1 
FZHD 0. 1 
CZLD 0. 1 
CZHD 0. 1 


Hall 

Mobility 

2 

(cm /V sec) 

208 

194 

212 

208 


Drift 
Mobility 
2 , 

(cm /V sec) 

226 

211 

230 

226 


Measured Hall mobilities for low-resistivity p-type 
silicon specimens. Drift mobilities were obtained 
assuming a Hall-to-drift mobility ratio of 0. 92 
(Reference 18). 
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3. 3 Resistivity Considerations 

During the preliminary stages of .this investigation, an apparent discrepancy 
was observed in some cases between four- point- probe resistivity measure- 
ments and resistivity determinations based on potential profiles. In partic- 
ular, differences between measurements on thick and thin samples were 
noted. In an attempt to resolve this discrepancy, careful four -point- probe 
and potential-profile measurements were made on four low- resistivity sam- 
ples, two of which were relatively thick (0. 32 cm) and the other two were 
relatively thin (~0. 03 cm). Results are shown in Table III and it is seen 
that the two methods yield resistivities that are in excellent agreement for 
thin specimens. For the thicker samples, four-point- probe values are a 
factor of 1. 07 larger than their potential- profile counterparts. Upon com- 
paring four-point-probe readings for thick and thin samples prepared from 
the same ingot, values for FZLD 0. 1 agree exactly. For CZHD 0. 1, the 
thin-sample value is a factor of 1. 07 larger than that for the thick specimen. 
In view of the fact that the discrepancies noted are rather small and that 
the two measurement techniques agree for thin samples, we decided to 
employ resistivities measured using a four-point probe for all of the thin 
bulk silicon samples under study. The source of observed differences is 
unknown at present, but determination of this source was considered un- 
important for the present investigation. 

TABLE III 

Resistivity 


Sample 

Designation 

Dimensions 
• cm 

F our - Point 
Probe 
Q-cm 

Potential Ratio of 

Profile Two Methods 
Q-cm 

CZHD 0. 1/689 

0. 60 x 1. 53 x 0. 32 

0. 106 

0. 099 

1. 07 

CZHD 0. 1/2 

0. 99 x 1. 99 x 0. 030 

0. 1 13 

0. 113' 

1. 0 

FZLD 0. 1/686 

0. 61 x 1. 54 x 0. 32 

0. 108 

0. 101 

1. 07 

FZLD 0. 1/10 

1. 0 x 2. 0 x 0. 031 

0. 108 

0. Ill 

0. 97 

Comparison of resistivity measurements made by four-point- 
probe and potential-profile methods for thick (0. 32 cm) and thin 
0. 03 cm) low-resistivity bulk silicon specimens. 



3. 4 Dislocation Density Determinations 

Dislocation densities were determined for the four low resistivity silicon 
ingots using the Sirtl etch method (ASTM Designation F-47-70). A chromic 
acid etch was employed and two slices from each ingot, designated A and B, 


17 



were examined. For a given ingot, the positions of slices A and B were 
chosen such that all other material to be studied (i. e. , material for prep- 
aration of bulk samples, solar cells, and slices for determination of 
oxygen concentration) was situated between them. This approach was 
employed to obtain an indication of variations in dislocation density over 
the region under investigation for a particular ingot. Results obtained 
are shown in Table IV. Measurements made by Dow Corning are shown 
for comparison. 


TABLE IV 

Dislocation Densities for Four Low Resistivity Silicon Ingots 


Ingot 

Designation 


Dow 
Corning 
Ingot No. 


Dislocation Density (cm ) 

Slice A Slice B Dow Corning 


FZHD 0. 1 
FZLD 0. 1 
CZHD 0. 1 
CZLD 0. 1 


S73-5277 
S71-1305 
S73-5278 
S 73-5276 


2. 5 x 10 4 

v.U 

<T 

3. 9 x 10 3 


2. 8 x 10 4 

vU 

O' 

0 

2. 9 x 10 3 


0 . 0 


1. 5 to 3 x 10 4 
0 

3 to 6 x 10 3 
0 


For FZLD 0. 1, a swirl pattern of etch pits was observed which apparently 
is not associated with the presence of dislocations. See text for discussion. 


Ingot FZLD was grown by Dow Corning to be dislocation free. However, we 
observed a swirl pattern of etch pits for slices A and B. These swirls cov- 
ered about 30% of the slice surface, with the remainder of the surface being 
free of etch pits. The local etch pit density in a swirl was ~1. 3 x 10^ cm" , 
and the average density (entire area) for both slices was ^4 x 10^ cm"^. 

Etch pit patterns such as those noted here have been observed in dislocation- 
free float-zone silicon by de Kock^’ ^ and Ciszek^ and have been attributed 
to vacancy clustering at a vacancy-oxygen complex. ^ Thus, the presence 
of etch pits does not necessarily indicate the presence of dislocations. Based 
on the dislocation density given by Dow Corning for this ingot (Table IV) and 
on the close correspondence of our observations of etch pits with those of 
others 19-2 1 on (ji s i 0 cation-free -material, we conclude that ingot FZLD 
most likely is free of dislocations. 
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3. 5 Vacuum Grease Experiments 

For the irradiations performed in this study, it was necessary to insure that 
bulk samples and solar cells were in good thermal contact with a hekt sink 
to prevent unwanted temperature rises during bombardment. As mentioned 
above, specimens were attached to a heat sink by using a small amount of 
Apiezon H vacuum grease. (This approach had previously been employed 
with success at JPL for solar cell irradiations). However, the question 
arose as to whether the grease could be completely, and readily, removed 
from a sandblasted surface. Additionally, we wished to determine whether 
this grease would affect the measured steady-state photoconductivity signal 
by changing the surface electrical and/or optical properties. We performed 
experiments which addressed these two questions. The following procedure 
was found to quickly, and apparently completely, remove Apiezon H: a sample 
with grease on it was dipped for a few seconds in a beaker of trichloroethylene 
contained in an ultrasonic cleaner; this treatment was followed by a few sec- 
onds in acetone, then methyl alcohol, and finally deionized water, all three 
of which were also in the ultrasonic cleaner. Examination under a micro- 
scope revealed no grease on the sample surface. To examine the effect of 
Apiezon H on photoconductivity signal, three 0. 3 mm thick samples were 
measured in the SS PC apparatus both before application of grease and after 
application of grease and subsequent cleaning of the surface using the above 
procedure. The average post-cleaning photoconductivity signal was found 
to be within a few per cent of the corresponding pre-grease measurement, 
and we concluded that Apiezon H can be employed to attach specimens to a 
heat sink without having a detrimental effect on subsequent lifetime measure- 
ments. 

3. 6 Generation Rate Calibration for the SSPC Apparatus 

Measurement of steady-state photoconductivity signal AV for a given spec- 
imen yields excess carrier density An through Equation 2. Lifetime is 
then obtained through Equation 3 assuming knowledge of the carrier gen- 
eration rate g. Absolute determination of diffusion lengths, and thus 
damage coefficients, depends on the accuracy to which g is known 
for the SSPC apparatus. A careful series of experiments were performed 
for this purpose, and the procedure followed and results obtained are 
briefly described in this section. 

Thirteen bulk silicon samples, listed in Table V, were employed in the 
calibration, all of which were prepared in the usual manner (see Section 
2.2). As a first step, photoconductivity decay techniques were used to 
determine low-level carrier lifetimes (designated T ^ in Table V). In- 
jection ratios at which PCD measurements were made ranged from ^10“^ 
for low-resistivity samples to ^5 x 10“^ for higher resistivity specimens. 

In most cases, particularly for CZ 2. 5 samples and FZLD 0. 1 samples, 
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T \ 


SAMPLE 

DESIGNATION 

RESISTIVITY 

ohm-cm 

T 

SS 

M-sec 

T pcd 

y.sec 

SS 

T 

pcd 

CZ 2. 5/9 

2. 58 

3. 59 

3. 59 

= 1. 0 

CZ 2. 5/10 

2. 54 

3. 64 

3. 58 

1. 02 

FZ 2. 5/1 

2. 51 

3. 82 

3. 60 

1. 06 

FZ 13/5 

14. 9 

3. 44 

3. 75 

0. 92 

FZ 13/6 

12. 4 

3. 77 

3. 75 

1. 01 

CZ 2. 5/1 

2. 52 

3. 71 

3. 67 

1. 01 

CZ 2. 5/2 

2. 54 

3. 90 

3. 66 

1. 07 

FZ 13/3 

12. 8 

3. 30 

3. 64 

0. 91 

FZ 13/4 

13. 2 

3. 34 

3. 54 

0. 94 

FZLD 0. 1/2 

0. 106 

6. 85 

6. 83 

1. 00 

FZLD 0. 1/3 

0. 104 

7. 36 

7. 16 

1. 03 

FZHD 0„ 1 /2 

0. 099 

5. 44 

5. 67. 

0. 96 

FZHD 0„ 1/3 

0. 099 

5. 72 

5. 95 

0. 96 


TABLE V 


Measured resistivities and effective minority carrier lifetimes for 
bulk silicon specimens. Resistivity measurements were made using 
a four-point probe and lifetimes were measured using both steady- 
state -photo conductivity (t ) and photoconductivity-decay ( T p C( j) 
techniques. 
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trapping effects were small. In a few cases (particularly for FZ 13), long- 
time-constant trapping effects were noted. Because lifetimes were much 
shorter than trapping times, a baseline shift occurred in the photoconductivity 
decay; recombination time constants were readily extracted from observed 
signals of this type. 

PCD lifetime measurements for CZ 2. 5 were judged to be the most accurate 
and we thus employed one of these samples (CZ 2. 5/9) as a standard for 
determining g. Steady- state measurements were then made in the SSPC 
apparatus at a constant chopped light intensity under a traps -filled condition 
(achieved with the aid of background lamps-- see Section 3. 7). A measure- 
ment of An was made for the standard sample, and generation rate was 
determined through the relation 


g = WT pcd . (13) 

The quantity An was also measured for the other samples using SSPC, and 
steady-state lifetime, T sg , then determined using the value of g obtained 
from the standard specimen. (Injection ratios due to chopped light in the 
steady- state measurements ranged from 1 to 3 x 10"^ for low resistivity sam- 

_ A 

pies and from 1 to 5 x 10 for higher resistivity specimens. ) The ratio of steady- 

state lifetime obtained in this fashion to t , is listed in Table V for each 

pcd 

sample and it is seen that all values lie crose to the ideal ratio of unity. 

With the exception of samples FZ 13/3 and FZ 13/5, all ratios are within 
7% of 1.0. Note that the quantity t /t j was determined for a resistivity 
range covering more than two orders of magnitude. The favorable results 
obtained suggest that our measurements of AV, V Q and p are quite accurate 
over the range studied and also -that the assumed mobility values are con- 
sistent with measured resistivities. 

The steady- state lifetime measurements listed in Table V were made at a 
generation rate (due to chopped light) of 1. 5 x 10^ pairs /cm^- sec, as 
were all other SSPC lifetime measurements in this investigation. With the 
present experimental set-up, the highest generation rate that can be achieved 
is a factor of ' v 2. 5 larger than this quantity. 

In calibrating the SSPC apparatus in terms of generation rate, it was neces- 
sary to exercise extreme care in the relative placement of all parts of the 
system (sample, filters, light source, etc. ). In particular, a considerable 
amount of time was spent in constructing a sample holder that would allow 
repeatable and reliable placement of a specimen with respect to the filters 
and the chopped light source. A slight rotational placement error, for 
example, results in an appreciable change in generation rate. Additionally, 
calibration of the SSPC apparatus was periodically checked with a set of 
unirradiated samples to ensure that g remained invariant during the course 
of the investigation. 
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It was considered desirable to check the calibration of the SSPC apparatus 
against an independent measurement method. Diffusion-length determinations 
by the steady- state surface photovoltage method were made on several sam- 
ples, and results of the SSPC--SPV comparison are discussed in Section 4.4. 

3. 7 Trapping Effects 

Measurements of relative lifetime versus injected-car rier density were made 
on representative bulk silicon samples with sandblasted surfaces upon receipt 
from Centralab (prior to the generation-rate calibration) to determine whether 
steady- state photoconductivity measurements were affected by minority-carrier 
trapping. The shape of lifetime- versus -exce^-density curves at low injection 
levels reveals whether trapping is important. 9 The SSPC apparatus was 
employed for these measurements in the following way. Excess density An, 
due to chopped-light illumination, was determined in the usual manner (Section 
2.1). We can express An as 


An = gT = g'lT , 


( 14 ) 


where g ' is a generation constant and I is illumination intensity. A signal 
proportional to I was measured using a silicon photodiode (refer to Figure 
1), and relative lifetime then determined by dividing An by this signal. 

Rather strong trapping effects were observed for all four low-resistivity 
ingots. Trapping was also noted for 2. 5 and 13 ohm- cm samples.* These 
results, particularly for 0. 1 ohm- cm material, were surprising to us and 
were not expected based on our previous observations. We have generally 
observed in the past that trapping effects become less important as resis- 
tivity decreases. Additionally, relatively thick Northrop-fabricated 2. 5 
ohm-cm samples with sandblasted surfaces had been observed earlier to 
exhibit relatively slight trapping. This observation suggested that the 
noted trapping effects were associated with surface properties. 

A few experiments were performed for the purpose of attaining insight into 
the nature of the observed trapping effects. Figure 4 shows measured rel- 
ative lifetime versus excess carrier density for sample FZHD 0. 1/1 under 
three conditions: a) original (as received from Centralab); b) after etching 
the surfaces with CP-4; c) after sandblasting the etched surfaces. Data 
are 'normalized to unity relative lifetime at an injected-carrier density of 
3 x 10^ cm~^ for convenience of comparison. It is seen that etching en- 
hanced trapping whereas a subsequent sandblasting diminished trapping 
effects compared to the original sample. The same experiment was per- 
formed for a 2. 5 ohm-cm sample and a similar result was obtained. There 
are at least two relevant statements that can be made based on these findings: 
a) Observed trapping effects are surface related since changing surface prop- 
erties changes trapping behavior; b) Results obtained following sandblasting 

#The observation of relatively strong trapping effects in SSPC measurements 
as compared to PCD can be explained plausibly in terms of significantly 
fewer filled traps in the transient measurement case. 
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Figure 4. Relative lifetime vs excess carrier density for a 0. 1 ohm- 
cm FZ Si sample after various surface treatments. (The 
’’original” case refers to a Centralab-fabricated sample. ) 
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(after etching) are consistent with findings on Northrop-fabricated sand- 
blasted samples in which trapping effects were observed to be considerably 
less important than for Centralab samples. Regarding the latter statement, 
an experiment was performed on sample CZLD 0. 1/1 in which lifetime 
versus excess density curves were obtained before and after a Northrop 
sandblasting (no etching). Results are shown in Figure 5 and it is seen 
that Northrop sandblasting did not alter the amount of trapping but merely 
decreased the amount of light absorbed in this Centralab-fabricated sample 
by about a factor of two due to reflectivity changes. 

The experiments described above aimed at determining the nature and origin 
of observed trapping effects are obviously incomplete. Further study was 
beyond the scope of the present work (particularly in view of the successful 
use of non-penetrating light to fill traps in Centralab-fabricated samples 
as described below). One might argue that all samples could be etched 
and then sandblasted to reduce trapping. However, it would be very dif- 
ficult to control such processing so that surface optical properties had 
the required sample -to- sample uniformity that is necessary in making 
accurate measurements of lifetime using the SSPC apparatus. It is quite 
likely that Centralab-fabricated samples have the required uniformity 
although they also have the trapping problem. We decided to employ un- 
modified Centralab bulk Si samples (in conjunction with non- penetrating 
light) to measure lifetimes and subsequently obtain damage coefficients. 

A common approach to separating trapping and recombination is to employ 
a background (unchopped) light to keep the traps filled which then permits 
observation of steady- state photoconductivity signals that are not influenced 
by trapping. As traps are filled by an increasingly intense background 
light, the observed signal will decrease until a saturation is reached which 
corresponds to the low-injection-level value of carrier lifetime. Further 
increases of background lamp intensity beyond this saturation point will 
not result in further changes in the photoconductivity signal. (Eventually, 
at excess densities higher than under consideration here, lifetime will 
become injection-level dependent and the above statement will no longer 
hold. ) Our first approach to filling traps was to employ penetrating back- 
ground illumination (i. e. , unchopped filtered light). However, it was 
found that at the highest available background lamp intensity all the traps 
were not yet filled. With an increasingly intense lamp, there is the problem 
of heating the Si filter. Bandgap changes with temperature will alter the 
spectrum of transmitted light, which is an undesirable effect. . Assuming 
that observed trapping was surface related, then the utilization of non- 
penetrating light to fill surface traps is appropriate. It was decided that 
non-penetrating light was probably the best solution to the trapping problem, 
and to implement this approach several small lamps were installed near 
the sample holder in the SSPC apparatus (see Figure 1). Light from 
these lamps does not pass through the Si filter and is thus highly absorbed 
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Figure 5. Relative lifetime vs. excess carrier density for a 0. 1 ohm- 
cm CZ Si sample before and after Northrop sandblasting. 

(The "original" case refers to a Centralab-fabricated sample. ) 
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near the surface of a bulk Si sample. The placement of the background 
lamps was such that uniform illumination of the front and rear surfaces 
of a given sample resulted. 

Steady- state photoconductivity measurements were then made for one un- 
irradiated sample of each type to determine whether the background lamps 
could fill traps completely. For all samples examined, a saturation con- 
dition was reached in which further increases in background (unchopped) 
lamp intensity did not cause a change in the chopped-light signal. This 
experiment demonstrated that nonpenetrating light can be used to separate 
trapping and recombination, thus permitting measurement of pre-irradiation 
lifetime for all specimens. 

In the course of examining background light effects, it was observed that 
the lamp intensities required to completely fill traps resulted in sample 
heating. Heating occurred rather slowly (~1° per min for the first few 
minutes), and we therefore adopted the approach of performing photo- 
conductivity measurements in a time ^ 1 min after the background lamps 
were turned on. Taking data in this manner presented no problem and 
allowed the undesirable effects of a sample temperature rise to be avoided. 

Upon discovering the unexpectedly strong trapping effects for bulk samples, 
it became apparent to us that even if trapping and recombination could be 
separated before irradiation using an appropriate background lamp, this 
would possibly be difficult to accomplish after irradiation. This situation 
would not arise, however, if the following reasoning holds. If it is as- 
sumed that the traps in question are located at the surface, that the sur- 
face recombination velocity is infinite, and that any increase in surface 
trap concentration due to irradiation is negligible compared to that already 
present, then the lamp intensity required to fill surface traps will be in- 
variant with irradiation. This will be true since the lifetime at the surface 
is assumed to be zero and thus that portion of the excess carrier density 
at the surface that is available for filling traps will not be affected by ir- 
radiation (surface lifetime cannot experience a further decrease). If the 
above model is not an accurate description of the physical situation for 
the samples in question, then it is possible that an increasingly intense 
background lamp would be required following irradiation to fill traps. 

This situation would occur if the fraction of injected carriers available 
for trapping at a given lamp intensity decreases as bulk lifetime decreases 
due to irradiation. In this case, one quickly reaches experimental 
limitations imposed by the choice of an appropriate background lamp. 

To assess our ability following irradiation to separate recombination and 
trapping using background lamps, a trial 1. 0-MeV electron irradiation 
was performed to a fluence of 3 x 10^ electrons /cm^, A careful post- 
irradiation examination of bulk silicon specimens revealed that with the 
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background-lamp set-up described above all traps could be filled, even 
though lifetimes were significantly degraded, thus permitting measure- 
ment of post-irradiation lifetime. Somewhat more intense illumination 
was required to fill traps compared to the pre-irradiation situation. It 
now appeared that we could separate trapping and recombination and ac- 
curately measure lifetimes in irradiated bulk samples. (As discussed in 
Section 4. 0, for certain cases in subsequent irradiations measured diffusion 
lengths appeared to saturate at high fluences, and this behavior is attributed 
to radiation-induced trapping effects. ) 

3. 8 Surface Recombination- -Experimental Considerations 

The employment of a high surface recombination velocity as a boundary 
condition for bulk silicon specimens was selected on the basis of repro- 
ducibility and stability of surface properties. Although techniques for 
producing low-s surfaces exist, it was felt that a considerable technical 
effort would have been necessary to ensure that a . selected technique would 
repeatably produce surfaces with a very low surface recombination velocity 
that did not vary with time. On the other hand, sandblasting is known to 
produce high-s surfaces and, as shown in Section 3. 1, s can vary from 
~10 4 cm/sec to infinity without significantly affecting measured photocon- 
ductivity signal. In this section, we describe a few preliminary experi- 
ments relating to surface recombination effects and also discuss consid- 
erations regarding the s-value employed for Centralab-fabricated bulk 
samples. 

3. 8. 1 Simulation of Zero and Infinite Surface Recombination Velocity 
with a Diffused Sample 

During the course of deciding on an appropriate boundary condition for bulk 
samples, an approach was suggested by NASA-Lewis personnel which, in 
principle, would allow one to achieve a condition of either s = 0 or s = ® 
on the same sample. Figure 6 shows a schematic representation of such a 
specimen. The diffused region between the ohmic contacts wraps around 
to the other side of the specimen and covers most of the non-contacted 
surface area. The sample is assumed to be prepared from p-type starting 
material with an n + diffusion. With lead A grounded, minority carriers 
created in the p-type material by illumination will be collected at the n + -p 
junction if they are within a diffusion length of it. Since the junction acts 
as a minority- carrier sink, a condition of s = °° is simulated. With lead 
A open-circuited, a condition of s = 0 is simulated. 

The above approach was experimentally examined using the SSPC apparatus. 
Centralab fabricated five n-on-p samples for us with the geometry shown 
in Figure 6. They utilized 10 ohm-cm p-type material for these specimens. 
Measurement of the steady- state photoconductivity signal was made with 
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Figure 6. Schematic representation of a diffused sample used 
to simulate conditions of s = 0 and s = ®>. 
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the diffused region open-circuited and with this region grounded for all 
five units. The following open-circuit-to-grounded photoconductivity- 
signal ratios were obtained for samples 1 through 5 respectively: 5. 4, 

3. 0, 1. 1, 3. 1, and 2. 1.' The larger open-circuit signal observed is con- 
sistent with the expectation of a low surface recombination velocity for 
such a configuration. 

Considering the result for sample #1, Figure 3 reveals that a value for 
t 6 ££/t of 0. 19 ( = 1/5.4) occurs for a bulk diffusion length of ~190 |im 
(using the s = 0 and s = ®> curves). For 10 ohm-cm p-type material, a 
pre-irradiation diffusion length of 190 p,m is reasonable for a solar cell, 
which suggests that the sample structure employed can be effectively 
utilized to simulate conditions of s = 0 and s = ®. However, it is obvious 
that additional study, both experimental and analytical, would be neces- 
sary before one could have confidence in such an approach to steady- state 
lifetime measurements on thin samples. In particular, a detailed analysis 
of the open-circuit condition is needed to demonstrate how closely a zero 
surface recombination velocity is approximated. We conclude that the 
diffused- sample approach to simulating conditions of s = 0 and s = 03 ap- 
pears promising, but additional work would be required to raise our level 
of understanding of this method to the point where it could be confidently 
employed in steady- state lifetime measurements. 

3. 8. 2 Measurement of the Thickness Dependence of Lifetime for a 
Sandblasted Bulk Sample 

An experiment was performed to examine the thickness dependence of 
measured minority-carrier lifetime in a bulk silicon sample with sand- 
blasted surfaces. The purpose of this study was to compare theory with 
experimental findings and thus evaluate the accuracy with which we could 
correct for surface recombination effects during steady-state-photocon- 
ductivity measurements of lifetime on solar-cell- size bulk silicon samples 
with sandblasted surfaces. A sample was prepared from ingot C-Z 2. 5 
(actual sample resistivity = 2.7 ohm-cm) and measured 30. 2 x 6. 89 x 
6. 82 mm after lapping to remove saw damage and subsequent sandblasting 
of all surfaces with 320 mesh silicon carbide. Ohmic contacts were pre- 
pared on the end faces using ultrasonic soldering techniques. Forward 
and reverse resistance measurements indicated no rectification at room 
and dry ice temperatures. Lifetime was measured using photoconductivity 
decay. The observed decay was exponential, yielding a value of 44 psec, 
and no trapping effects were observed. Such a situation is ideal for tran- 
sient measurements of lifetime. Steady-state photoconductivity measure- 
ments were also made, as discussed below. 

Figure 7 shows measured lifetime (PCD) versus sample thickness. The 
original sample was cut into two pieces which were 4. 1 and 2. 4 mm thick 
after lapping to remove saw damage and then sandblasting. Lifetime for 
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(oasri) 3WIJL3Jn 

Figure 7. Measured lifetime (PCD) versus sample thickness for a 2. 7 
ohm-cm boron-doped silicon sample. Shown for comparison 
are analytical results for two values of surface recombination 
velocity using a one-dimensional model. 
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sample A (4. 1 mm) was measured on two consecutive days to check re- 
peatability. A slight difference was observed (~6%) which probably re- 
flects experimental error or possibly a small temperature difference. 

To obtain samples C and D, sample B was cut into two equal- thickness 
pieces. As shown in Figure 7, measured lifetimes for these samples 
are in reasonable agreement. To obtain sample E, sample C was lapped 
down to half its original thickness and then sandblasted. Sample F was 
obtained by lapping sample E. Following lapping and a light sandblasting 
of the surfaces, the sample thickness was 0. 31 mm, which is nearly the 
same as the thickness of bulk samples prepared by Centralab and utilized 
here for damage coefficient determinations. To examine for possible 
changes in effective lifetime due to changes in surface properties with 
time, lifetime measurements were repeated five weeks later for several 
samples (A, D, and E). Within a few per cent, lifetime was found to be 
invariant over this time period. 

For purposes of analysis, Equation 6 was employed. Because this ex- 
pression is one-dimensional (i. e. , assumes two of the sample dimensions 
are large compared to the third), it does not apply for large sample thick- 
nesses, and this point should be kept in mind when comparing experimental 
anc^ analytical results (Figure 7). A value for electron mobility of 1150 
cm /V-sec was assumed, yielding a D value of 29. 8 cm^/sec. For bulk 
lifetime t, the value measured for the original sample (44 p-sec) was em- 
ployed, resulting in a bulk diffusion length of 362 |im. Effective lifetime 
was calculated as a function of sample thickness for two values of surface 
recombination velocity, and results are shown in Figure 7. For s = 00 , 
agreement with experiment is good for 0. 31, 0. 5, 1, and 2. 4 mm specimens. 
For the thicker samples, Equation 6 is not applicable. Note that the mea- 
sured lifetime at 0. 31 mm is only 6% of that for the original sample, and 
for s = ® the theory predicts a value for effective lifetime that is in ex- 
cellent agreement. For s = 10^ cm/sec, however, the predicted value at 
0. 31 mm is a factor of ^5 too large. The assumed bulk value of 44 iisec 
could possibly be too small due to surface effects even on the largest 
sample. However, a larger value of t would result in an even greater 
discrepancy between theory and experiment for s = 10^ cm/sec. Con- 
sideration of experimental results and analysis leads us to conclude that 
the surface recombination velocity for a sandblasted surface is indeed 
large, as expected, and that measurements agree closely enough with 
theory to allow a reasonably accurate correction for surface recombination 
effects in the present investigation. 

Steady- state photoconductivity measurements were also made versus 
sample thickness using the SSPC apparatus. For these measurements, 
the thickness of the Si filter was 2. 5 mm. Since this filter thickness is 
< the sample thickness for the original sample and for samples A and B, 
a nonuniform excess carrier density results in these samples, which 
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complicates data interpretation. Therefore, measurements on relatively 
thin samples (C through F) are more readily interpretable. Upon compar- 
ing [relative t (1 mm) /relative t (0. 5 mm)] for steady-state measurements 
with the corresponding absolute quantity for transient data (Figure 7), rea- 
sonable agreement (within ~15%) was found. Once again, the indication is 
that the experimental situation is understood well enough to be able to apply 
the theory directly in accounting for the effects of surface recombination. 

3. 8. 3 Surface Recombination Velocity for Bulk Specimens 

As discussed in Section 3. 1, Equation 10 was employed to correct measured 
effective lifetimes for the effects of surface recombination. An iterative 
procedure was utilized with this equation to obtain actual bulk diffusion 
length. As a starting point, we assumed s = co and then substituted various 
values of L into Equation 10. This procedure did not work for any of the 
Centralab -fabricated bulk samples; no finite value of L would yield a value 
f°r T e g as large as those measured, assuming infinite s. This result sug- 
gested that s was actually less than infinite for Centralab -fabricated sam- 
ples. As discussed above, measurements of lifetime versus sample thick- 
ness for a CZ 2. 5 specimen with Northrop-sandblasted surfaces were made, 
and excellent agreement was found between analysis and experiment when 
infinite s was assumed. This finding strongly indicated that the Northrop 
surface treatment results in an infinite surface recombination velocity. 
Centralab also sandblasted sample surfaces after preparation, but subse- 
quently it became apparent that there were differences between the surface 
optical properties of Northrop-sandblasted samples and Centralab samples. 
At a given light intensity, measured excess carrier density was significantly 
larger in the Centralab case, indicating more energy absorption in the sam- 
ple bulk. This observation is also consistent with less recombination of 
excess carriers at the surface for Centralab samples (fewer surface re- 
combination centers) and thus a smaller s. More concrete evidence for a 
smaller s was sought, along with an appropriate method for determining s 
for Centralab samples. The procedure followed is now described. 

Samples from three higher-resistivity Northrop- supplied ingots were also 
fabricated by Centralab (FZ 2. 5, CZ 2. 5, and FZ 13). For these ingots, 
we had on hand accurate (negligible trapping effects) PCD lifetime mea- 
surements for large bulk samples where surface effects can be neglected. 
For FZ 2. 5 and FZ 13, we had fabricated a large number of bulk samples 
on another program three years ago and then measured pre -irradiation 
values of t. These existing data were used to obtain an average bulk life- 
time for each of the two ingots under consideration. Next, all measured 
pre-irradiation values of t ££ for each of these ingots were averaged, as 
were sample thicknesses, resistivities, and assumed mobilities. This 
procedure thus yielded. average values for the following quantities for each 

of two ingots based on a rather large number of samples: t , r (and 

eff 
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hence L), T, and D. Next, these average values were employed in Equation 
10 and the value of s that would then satisfy the equation was then determined. 
For FZ 2. 5 and FZ 13, we obtained s = 8. 3 x 10 3 cm/s,ec and s = 8. 0 x 10 3 
cm/sec, respectively. The same procedure was followed for ingot CZ 2. 5 
except that only one bulk lifetime value was available ( t = 44 jjisec). Using 
average values for T e fp T, and D along with this value of t, a value of 
s = 8. 5 x 10 3 cm/sec resulted. The evidence thus strongly indicated that 
the surface recombination velocity for Centralab-fabricated samples is 
slightly less than 10^ cm/sec. Referring to Figure 3, it is seen that the 
difference between analytical results for s = 10 4 and s = 00 is relatively 
small. However, in attempting to correct values of . t ^ for the effects of 
surface recombination in the present investigation, one cannot employ 
s = oo for Centralab samples. Based on the above results, a value of s = 

“3 

8. 3 x 10° cm/sec was assumed for all bulk silicon samples under study, 
and t was determined using Equation 10 in the manner described above. 

As discussed above, in terms of obtaining it is not necessary to per- 
form a correction for surface effects. This correction is obviously neces- 
sary, however, if one wishes to determine actual pre-irradiation lifetimes 
and diffusion lengths. 

3. 9 Determination of Oxygen Concentration 

An attempt was made to determine oxygen concentration in the four 0. 1 
ohm-cm silicon ingots employed in this investigation by using the infrared 
absorption method (ASTM Designation F-121-70 T). To perform this mea- 
surement for FZ material, a specimen with negligible oxygen concentration 
is required as a reference. A piece of oxygen-free material was supplied 
to us by Dow Corning and this piece, along with four low-resistivity slices 
(one per ingot) were mechanically polished according to specified ASTM 
procedures. For comparison, a slice from ingot CZ 2. 5 was also pre- 
pared for determination of oxygen content. Initial measurements were 
made using Beckman IR-4 and IR-20A spectrophotometers. Subsequently, 
a Beckman Model 42 50 IR spectrophotometer was employed. 

Measurements at room temperature on CZ 2. 5 revealed the expected ab- 
sorption peak at 9 (rm due to oxygen. However, in room temperature mea- 
surements on all four low- resistivity slices no transmission was observable 
for wavelengths ^5 um, thus preventing the observation of 9-ji*m absorption. 
Comparison of these negative results with the findings for CZ 2. 5 suggests 
that the problem with measurements on 0. 1 ohm-cm material is lack of 
transmission due to free carrier absorption. According to ASTM Desig- 
nation F-121-70 T, free carrier absorption should not be a problem for 
resistivities greater than 0. 01 ohm-cm in the method being employed. 

We also made measurements on low-resistivity Czochralski material at 
liquid nitrogen temperature for the purpose of freezing out majority car- 
riers and thus decreasing free carrier absorption. Per cent transmission 


33 



increased at wavelengths considerably less than 9 (J,m, but once again no 
transmission was observed at 9 p,m. Other workers have also been unable 
to observe transmission through 0. 1 ohm-cm silicon slices and were thus 
unsuccessful in determining oxygen content. ^ 

2 5 

It was subsequently learned that the ASTM statement regarding applica- 
bility of the method for resistivities 5:0. 01 ohm-cm is in error and should 
actually be ^0. 1 ohm-cm. Measurements on 0. 1 ohm-cm Czochralski 
material are quite difficult and apparently are not possible on float-zone 
material of that resistivity. To make measurements on Czochralski ma- 
terial, one apparently must violate ASTM procedures and make the test 
specimen thin enough (^0. 5 mm) to permit transmission. Additionally, 
the slit on the spectrophotometer should be opened as wide as permissible, 
given the 32 cm -1 room temperature line width, scale expansion should be 
employed, and a slow scan should be performed. 

The above approach was employed to measure the oxygen concentration in 
a thin (0. 35 mm) polished CZHD 0. 1 slice using a Beckman Model 4250 
dual-beam spectrophotometer. Transmission was observed at 9 Lim, and 
the oxygen absorption was a small perturbation on the free carrier absorption 
of about 50 cm” 1 . The observed absorption coefficient of 0. 9 cm” 1 , cor- 
rected for lattice band absorption, corresponds to an oxygen concentration 
of about 4 x 10^ cm“^. The observed linewidth exceeded the expected 
32 cm” 1 by 40%, however, so that the observed peak absorption is an under- 
estimate in view of the excessive slitwidth having been employed. 

3. 10 Pre -irradiation Properties of Bulk Samples and Solar Cells 

Table VI lists bulk silicon samples employed in the present investigation. 
Shown are resistivities and pre-irradiation values for effective lifetime, 
actual lifetime, and actual diffusion length. Also listed are irradiating 
particle energies for each sample and additional information regarding 
special conditions for certain specimens. As indicated, post-irradiation 
annealing studies were performed for selected samples at either room 
temperature (RT) or 60°C. Additionally, two samples were irradiated 
in the dark with 10-MeV protons to examine the effect of illumination on 



Table VII lists pre-irradiation properties of solar cells obtained using a 
solar simulator. Cells with the designation M dose monitor” were irradiated 
for the purpose of obtaining a rough check on fluences received by comparing 
cell degradation with that expected based on previous findings (such as Ref- 
erence 11). For the cells listed in Table VII, solar- simulator current- 
voltage characteristics were obtained following each bombardment, and 
data are presented in Section 4. 0. Also in that section are results of 
diffusion length measurements (using a Co-60 source) for these cells and 
other units. 
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Pre- Irradiation 


Sample 

Resistivity 

T 

eff 

,T 

Designation 

ohm -cm 

psec 

psec 

FZLD 0. 1/1 

0. 106 

7. 21 

37. 1 

4 

0. 110 

7. 09 

25. 5 

5 

0. 105 

7. 55 

46. 5 

6 

0. 105 

7. 54 

37. 1 

8 

0. 108 

7. 22 

27. 0 

10 

0. 108 

7. 08 

34. 2 

11 

0. 109 

7. 26 

2 7. 5 

12 

0. 108 

7. 49 

34.2 

13 

0. 105 

7. 60 

52. 9 

14 

0. 107 

7. 43 

36. 2 

15 

0. 109 

6. 71 

22. 7 

16 

O’. 105 

7. 25 

40. 6 

17 

0. 105 

7. 12 

30. 5 

18 

0. 107 

7. 38 

41. 5 

19 

o. no 

6. 90 

24. 0 

20 

0. 105 

7. 30 

42. 4 

21 

0. 107 

6. 96 

33. 5 

22 

0. 105 

7. 56 

54. 3 

FZHD 0. 1/4 

0. 100 

5. 69 

15. 3 

5 

0. 095 

6. 03 

19. 1 

6 

0. 097 

5. 97 

20. 0 

7 

0. 096 

5. 93 

20. 4 

8 

0. 099 

5. 87 

18. 7 

11 

0. 099 

5. 99 

20. 7 

12 

0. 096 

5. 93 

20. 4 

13 

0. 093 

6. 37 

25. 9 

14 

0. 096 

6. 14 

21. 9 

15 

0. 097 

5. 73 

17. 7 

16 

0. 099 

5. 92 

19. 6 

17 

0. 096 

5. 85 

18. 5 

18 

0. 095 

5. 99 

21. 0 

19 

0. 096 

6. 24 

23. 0 

20 

0. 094 

6. 02 

21.4 

21 

0. 095 

6. 03 

20. 7 

22 

0. 095 

6. 17 

22. 2 

23 

0. 094 

6. 10 

22. 3 

24 

0. 094 

6. 08 

21.2 

25 

0. 097 

5. 73 

18. 3 


TABLE VI 


L 

Irradiating 

Particle 

Comments / 
Special 

pm 

Energy, MeV 

Condition 

201 

0. 5 

RT anneal 

167 

1 1 

- 

225 

1 1 

- 

201 

1 1 

- 

172 

1. 0 

RT anneal 

193 

1 1 

60°C anneal 

173 

1 1 

- 

1-93 

n 

- 

240 

1 1 

- 

199 

2. 5 

RT anneal 

157 

1 1 

- 

210 

1 1 

- 

182 

1 1 

- 

213 

10 

RT anneal 

162 

1 1 

60°C anneal 

215 

1 ! 

- 

191 

1 1 

- 

243 

! 1 

- 

129 

10 


146 

0. 5 

RT anneal 

148 

l 1 

- 

149 

I 1 

- 

143 

t t 

- 

149 

1. 0 

RT anneal 

149 

1 1 

60°C anneal 

168 

10 

unilluminated 

155 

1. 0 


139 

1 1 

- 

146 

1 1 

- 

142 

2. 5 

RT anneal 

151 

1 1 

- 

158 

1 1 

- 

153 

} 1 

- 

150 

10 

RT anneal 

156 

n 

60°C anneal 

156 

1 1 

- 

152 

f f 

- 

141 

1 1 



Pre-irradiation lifetimes and diffusion lengths for bulk silicon specimens of 
various resistivities. 
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Pre -Irradiation 


Sample 

Resistivity 

T eff 

T 

L 

Irradiating 

Particle 

Comments / 
Special 

Designation 

ohm- cm 

usee 

ixsec 

p,m 

Energy, MeV 

Condition 

CZLD 0. 1/3 

0. 117 

1. 45 

1. 84 

44. 7 

0. 5- 

RT anneal 

5 

0. 112 

1. 62 

2. 15 

48. 4 

It 

- 

6 

0. Ill 

1. 46 

1. 89 

45. 4 

It 

- 

7 

0. 113 

1. 52 

1. 98 

46. 5 

1 1 

- 

8 

0. 110 

2. 00 

2. 76 

54. 8 

10 

- 

9 

0. 112 

1. 66 

2. 16 

48. 5 

1. 0 

RT anneal 

10 

0. Ill 

1.. 74 

2. 30 

50. 0 

' II 

60°C anneal 

11 

0. 114 

1. 76 

2. 34 

50. 5 

1 I 

- 

12 

0. 116 

1. 60 

2. 10 

47. 8 

II 

- 

13 

0. 114 

1. 78 

2. 36 

50. 7 

1 1 

- 

14 

0. 110 

1. 86 

2. 54 

52. 6 

2. 5 

RT anneal 

15 

0. 113 

1. 65 

2. 17 

48.6 

1 1 

- 

16 

0. 115 

1. 68 

2. 19 

48.8 

II 

- 

17 

0. 115 

1. 82 

2. 42 

51.3 

II 

- 

18 

0. 114 

1. 83 

2. 48 

52.0. 

10 

RT anneal 

19 

0. 112 

1. 66 

2. 21 

49. 0 

II 

60°C anneal 

22 

0. 113 

1. 73 

2. 30 

50. 0 

1 1 

- 

23 

0. 112 

1. 68 

2.21 

49. 0 

1 1 

- 

24 

0. Ill 

1. 51 

1. 96 

46.2 

II 

- 

CZHD 0. 1/1 

0.112 

2. 38 

3. 53 

62.0 

0. 5 

RT anneal 

5 

0. 108 

2. 33 

3. 59 

62. 5 

1 1 

- 

6 

0. Ill 

2. 20 

3. 25 

59. 5 

1 1 

- 

7 

0. 109 

2. 04 

2. 98 

57. 0 

1 1 

- 

8 

0. 108 

2. 45 

3. 85 

64.8 

1. 0 

RT anneal 

9 

0. 105 

2. 75 

4. 74 

71. 8 

I 1 

60°C anneal 

10 

0. 109 

2. 18 

3. 26 

59. 6 

I I 

- 

11 

0. 103 

2. 49 

4. 00 

66.0 

II 

- 

12 

0. 108 

2. 74 

4. 63 

71. 0 

1 I 

- 

13 

0. 109 

2. 38 

3. 63 

62. 8 

2. 5 

RT anneal 

14 

0. 109 

2. 01 

2. 90 

56. 2 

1 I 

- 

15 

0. 112 

2. 04 

2. 93 

56. 5 

1 I 

- 

16 

0. 112 

2. 08 

3. 00 

57. 2 

1 1 

- 

17 

0. Ill 

2. 00 

2. 88 

56. 0 

10 

RT anneal 

18 

0. 112 

2. 04 

2. 90 

56. 2 

1 1 

60°C anneal 

19 

0. 109 

2. 45 

3. 74 

63. 8 

1 I 

- 

20 

0. 108 

2. 48 

4. 06 

66. 5 

1 1 

- 

21 

0. 112 

2. 02 

2. 90 

56. 3 

II 

- 


TABLE VI (continued) 


36 



Pre- Irradiation 




T 

T 

L 

Irradiating 

Sample 

Resistivity 

eff 



Particle 

Designation 

ohm- cm 

usee 

usee 

|im 

Energy, MeV 

FZ 2. 5/2 

2. 57 

3. 63 

25. 1 

288 

0. 5 

4 

2. 56 

3. 70 

42. 6 

375 

1 1 

5 

2. 43 

3. 95 

73. 2 

485 

1 . 0 

6 

2. 34 

4. 16 

140 . 

670 

1 1 

7 

2.47 

3. 77 

38. 2 

353 

2. 5 

8 

2. 50 

3. 97 

72. 1 

485 

? 1 

9 

2. 60 

3. 82 

46. 1 

390 

10 

10 

2. 43 

3. 55 

48. 6 

407 

1 1 

11 

2. 46 

4. 01 

87. 7 

535 

1 1 

CZ 2. 5/3 

2. 57 

3. 77 

24. 6 

285 

0. 5 

4 

2. 54 

3. 91 

30. 4 

315 

1 1 

5 

2. 57 

3. 78 

26. 3 

295 

1 . 0 

7 

2. 60 

3. 95 

37. 1 

350 

1 1 

8 

2. 53 

3. 76 

24. 5 

283 

2. 5 

11 

2. 39 

4. 30 

40. 4 

360 

10 

12 

2. 38 

4. 55 

63. 1 

450 

1 t 

13 

2. 41 

4. 16 

116 

610 

2. 5 

15 

2. 45 

4. 00 

86. 1 

530 

10 

16 

2. 53 

4. 04 

99. 6 

570 

1 1 

CZ 6/1 

6.49 

3. 98 

11.4 

198 

0. 5 

2 

5. 57 

4. 30 

14. 4 

226 

1 . 0 

3 

5. 90 

4. 35 

15. 5 

2 34 

2. 5 

4 

5. 93 

3. 29 

7. 27 

157 

10 

5 

6. 63 

3. 69 

9. 32 

177 

It 

6 

5. 93 

4. 46 

16. 9 

239 

1 1 

FZ 13/2 

13. 5 

3. 65 

48. 6 

411 

0. 5 

7 

15. 4 

3. 92 

117 

638 

1 1 

8 

13. 7 

3. 47 

29. 5 

320 

1 . 0 

9 

12. 6 

3. 76 

56.8 

440 

II 

10 

12. 5 

3. 55 

48. 6 

407 

10 

11 

13. 6 

3. 64 

87. 3 

551 

1 1 

13 

13. 4 

3. 69 

131 

675 

2. 5 

14 

12. 6 

3. 78 

126 

655 

1 1 

15 

12. 7 

3. 76 

55. 5 . 

435 

10 

18 

12. 3 

3. 82 

119 

630 

It 


TABLE VI (continued) 


Comments / 

Special 

Condition 


unilluminate d 
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Pre- Irradiation 


Solar Cell 

I 

sc 

V 

oc 

P 

max 

Irradiating 

Particle 


Designation 

mA 

mV 

mW 

Energy, MeV 

Comments 

FZLD 0. 1/SC2 

65. 0 

617 

28. 8 

1. 0 

1 cmx 2 cm 

4 

64. 5 

615 

28. 0 

10 

it 

FZHD 0. 1/SCI 

61. 2 

612 

. 27. 5 

1 . 0 * 

M 

2 

62. 5 ' 

613 

28. 4 

10 

1 1 

CZLD 0. 1/SC 2 

53. 7 

605 

23. 3 

1. 0 

H 

3 

54. 2 

600 

22. 0 

10 

1 1 

CZHD 0. 1/SC2 

58. 5 

605 

24. 9 

1.0 /;• 

1 1 

3 

58. 2 

605 

25. 5 

10 

1 1 

CZ 2. 5/SC2 

65. 4 

578 

27. 6 

1. 0 

1 1 

3 

66. 5 

578 

28. 1 

10 

1 1 

CZ 9/SC2 

144 

548 

57. 2 

0. 5 

2 cmx 2 cm 
(dose monitor) 

3 

145 

550 

57. 5 

1. 0 

it 

4 

145 

550 

57. 2 

2. 5 

1 1 

CZ 2. 5/SC4 

67. 2 

578 

28. 5 

10 

1 cmx 2 cm) 
(dose monitor) 


TABLE VII 

Pre-irradiation electrical properties of n-on-p silicon solar cells of various 
base resistivities (AMO, 140 mW/cm^, ' yj 27°C). 
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SECTION 4. 0 


RESULTS 


In this section, experimental results of this investigation are presented. 
Included are data showing the degradation of diffusion length for electron- 
and proton-irradiated bulk samples and solar cells of various resistivities. 
Damage coefficients are determined as a function of resistivity and device 
and bulk behavior are compared. Results of annealing studies on irradiated 
bulk specimens are given. Degradation with irradiation of solar cell param- 
eters obtained using a solar simulator is presented. A comparison is made 
of data obtained using both steady- state photoconductivity and steady- state 
surface photovoltage methods. 

4. 1 Results of Electron Irradiations 

Tables VIII through X list diffusion length as a function of fluence for bulk 
silicon samples irradiated with 0. 5-, 1. 0-, and 2. 5-MeV electrons, re- 
spectively. For the purpose of determining damage coefficients the quantity 
(L“^ - Lq"^) was plotted versus fluence, and a unity slope fit to data on a 
log-log plot yields a value for (refer to Equation 1). 

4. 1. 1 0. 5-MeV Electrons 

Figure 8 presents data for 0. 5-MeV electron-irradiated specimens of various 
resistivities. With the exception of data for CZ 6/1 and CZ 9/SC2, all data 
points shown are averages for samples of a given type (refer to Table^ VIII). 
Unity-slope lines fit most of the data, with a few exceptions. Data for 
CZLD 0. 1 and CZHD 0. 1 at the lowest fluence appear to be inaccurate, which 
is probably attributable to experimental error since only a very slight de- 
gradation in L occurred at this fluence. At the highest fluences, a saturation 
is beginning to occur for CZLD 0. 1 and CZHD 0. 1. (This effect was more 
markedly observed at higher energies, and these findings are discussed 
below. ) Data for FZ 13 cannot be fit with a line of unity slope. The reason 
for this is unknown, but is possibly attributable to experimental errors as- 
sociated with the relatively small amount of damage introduced or to an error 
in measured L~. To obtain for FZ 13 samples, values of (L“^ _ Lq - ^) 
at the highest fluence were employed since these values are the least affected 
by errors in Lq. For other resistivities, corresponds to the dashed- 
line fits in Figure 8. 


Figure 9 presents K-^ vs resistivity for 0. 5-MeV electron-irradiated speci- 
mens. It is seen that the bulk data fall into two groups and both can be fit 
with a line of slope equal to -0. 67. The fit for CZLD 0. 1, CZHD 0. 1, and 
CZ 6 yields values of K-^ that are larger by a factor of ~ 3 than those obtained 
from the fit for FZLD 0. 1, FZHD 0. 1, FZ 2. 5, CZ 2. 5, and FZ 13 if com- 
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Sample Pre-irradiation 3 x 10 e/cm 
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Diffusion length as a function of fLuence for 0. 5-MeV electron-irradiated bulk silicon specimens. 
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TABLE IX 

Diffusion length as a function of fluence for 1. 0-MeV electron-irradiated bulk silicon specimens. 


Sample Pre -irradiation 3 x 10 e/cm 
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Diffusion length as a function of fluence for 2. 5-MeV electron-irradiated bulk silicon specimens. 



ELECTRONS 



(,UI0) 0-1-1 


-2 -2 

The quantity (L - ) vs fluence for 0. 5-MeV electron- 

irradiated bulk silicon specimens. A data point obtained 
from an irradiated solar cell is also shown. 
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FLUENCE (e/cm 






O FZLD 0.1 
□ FZHD 0.1 
A CZL0 0.1 
0 CZHDO.l 



RESISTIVITY (ohm-cm) 




parison is made at a given resistivity. Except for the CZ 2. 5 result, these 
findings are consistent with observations of others^' ^ that FZ solar cells 
are more tolerant of electron irradiation than are CZ units. (This observa- 
tion and others presented below for different energies lead us to suspect that 
ingot CZ 2.5 was actually FZ material. However, infrared absorption measure- 
ments on a sample from this ingot revealed an absorption peak at 9 M>m, which 
is characteristic of a large oxygen content. Thus, ingot CZ 2.5 appears to 
be CZ material. It is not clear at present why damage coefficients for 
specimens from this ingot exhibited behavior consistent with that for FZ 
material. ) A larger K-^ for CZ material indicates that the presence of oxygen 
enhances production of those defects primarily responsible for diffusion- length 
degradation. 

Figure 9 reveals that Kj^- values for FZLD 0. 1 and FZHD 0. 1 are nearly 
identical. Damage coefficient values for CZLD 0. 1 and CZHD 0. l.are also 
in good agreement with each other. This observation indicates that there is 
no significant dependence of on dislocation density for 0. 1 ohm-cm 
material bombarded with 0. 5-MeV electrons. (Additional discussion relating 
to this point is given in the Appendix. ) 

One solar cell data point is shown in Figure 9. It is seen that the Kj^ value 
for CZ 9/SC2 is a factor of ~ 4 larger than that expected on the basis of 
findings for CZ bulk samples. Further comparisons of solar-cell and bulk 
findings are discussed below. 

As discussed in Section 2. 2, an incident 0. 5-MeV electron degrades in energy 
to ~ 0. 44 MeV by the time it reaches the center of a bulk sample with a thick- 
ness of 0. 3 mm. Thus, the bulk damage effects presented above are more 
appropriately associated with an average electron energy of ~ 0. 44 MeV 
rather than the incident particle energy. (This consideration can account for a 
portion of the difference between solar cell and bulk material damage coeffi- 
cients (at 0. 5 MeV) because the SSPC measurement technique samples damage 
throughout the entire test specimen whereas solar cell measurements are 
presumably only affected by damage within a diffusion length of the n~^p junction). 

4. 1. 2 1. 0 -MeV Electrons 

Figure 10 presents data for 1.0-MeV electron-irradiated bulk samples and solar 
cells. Low-resistivity specimens exhibit linear degradation with fluence at low 
fluences, but exhibit saturation at high fluences (i. e. , an apparent sharply 
decreased defect introduction rate). We defer additional discussion of ob- 
served saturation effects to Section 4. 3. Higher resistivity samples also 
indicate a deviation from linearity at high fluences. (At the highest fluence, 
the onset of carrier removal was noted for the higher resistivity specimens. 

The resistivity increase in FZ 13 samples was ~ 7% and was ~ 3% for FZ 2.5 
and CZ 2. 5. ) Data for several solar cells (FZLD 0. 1/SCI and SC2, CZLD 
0. 1/SCI and SC2) were obtained at two fluences and the observed degradation 
with fluence is close to linear for these two cell types. However, the two 
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The quantity (L - Lq ) vs. fluence for 1. 0-MeV electron- 
irradiated bulk silicon specimens and silicon solar cells. 
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solar-cell data points for ingot CZ 2. 5 (SCI and SC2) cannot be fit with a unity 
slope line. The degradation rate is apparently lower at the higher fluence 
examined. This result is qualitatively consistent with the observations of 
Meulenberg and Treble^® who noted a fluence dependence of for 1. 0-MeV 
electron-irradiated solar cells (~ 9 ohm-cm). 

Data corresponding to linear degradation regions (i. e. , good fits to unity- 
slope lines in Figure 10) were used to determine damage coefficients for bulk 
samples and Figure 11 presents the results along with comparison data (for 
solar cells) obtained by other workers'. 28-30 p or the present solar-cell 
findings, data in Figure 11 represent K^-values obtained at a single fluence. 

A line with a slope of -0. 63 fits several groups of data points well. Findings 
for FZ bulk samples, and also CZ 2. 5, are fit well by such a line. Bulk 
data for CZLD 0. 1, CZHD 0. 1, and CZ 6 are fit well by a parallel line. 

FZ material is again observed to be more radiation resistant than CZ, with 
the bulk CZ-to-FZ K ^ ratio being ~ 1 . 9, as compared to ~ 3 for 0. 5-MeV 
electron irradiation. 

Two data points taken from solar cell findings of Meulenberg and Treble^® 
are shown in Figure 11. They observed a fluence-dependent value, and 
both high- and low-fluence values are given in the figure. It is seen that this 
range of values overlaps the present solar cell data and bulk sample data. 

We further note that the fluence dependence of K-^ presently observed for 
CZ 2. 5 cells (SCI and SC2) encompasses a variation of by a factor of 
~ 1. 7, whereas the variation of for data of Meulenberg and Treble is a 
factor of ~ 3. This difference is presumably attributable to the fact that their 
observations were made over a wider fluence range than that examined here. 

Also shown in Figure 11 are four data points obtained from the work of 
Downing et al. ^ The dashed-line fit shown through their data also fits a 
number of the solar cell data points obtained in the present investigation. 

Data of Faith^® are also presented, and the resistivity dependence of Kj^ is 
seen to be about the same as for other data in Figure 11, although his K^- 
values are somewhat larger. 

Low- resistivity solar cell data in Figure 11 are observed to be consistent with 
their bulk- sample counterparts in that values of K ? for CZ cells are a factor 
of ~ 1. 9 larger than those for FZ units. Additionally, as was the case for 
0. 5-MeV electron-irradiated specimens, there appears to be no significant 
dependence of on dislocation density for low -resistivity bulk samples 
or solar cells. We further note that K^-values for low- resistivity solar cells 
are a factor of ~ 1. 8 larger than their bulk-material counterparts. For the 
higher- resistivity cases, solar -cell -to -bulk -mate rial ratios range from 

~ 2 to ~ 3. 2. 

4. 1. 3 2. 5-MeV Electrons 

Data for 2. 5-MeV electron-irradiated specimens are shown in Figure 12, 
and saturation effects are observed at all resistivities. This phenomenon is 
particularly strong for low-resistivity specimens, and to determine for 
such samples we employed the values of (L - ^ - Lq"^) obtained at the lowest 
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- 2 ? 

The quantity (L, - Lq~^) vs. fluence for 2. 5-MeV electron-irradiated 

bulk silicon specimens, A data point obtained from an irradiated solar 
cell is also shown. 


Figure 12. 




13 2 

fluence ( 3 x 10 e/cm ). Damage coefficients obtained in this manner can 
be considered lower limits to actual values. As was the case for 1.0-MeV 
electron irradiation, carrier removal was observed here also. At the highest 
fluence, resistivity increases of ~'25% for FZ 13 and ~ 6 % for FZ 2.5 and 
CZ 2. 5 were noted. 

Figure 13 presents damage coefficient vs resistivity and it is seen that most 
of the data are fit well by a line with a slope of -0. 67. Based on the CZ 6/3 
data point, a CZ-to-FZ Kj^ ratio of ~ 1.6 is obtained. Low-resistivity samples 
do not exhibit such a ratio, but those data are most likely influenced somewhat 
by the saturation effect. As observed in Figures 8 and 10 , differences between 
FZ and CZ material diminish as saturation is approached. We further note that 
the one solar cell data point in Figure 13 (CZ 9/SC4) would yield a solar-cell - 
to-bulk-material K-^ ratio of ~ 1.4 if one were to extrapolate a line of -0. 67 
slop through the CZ 6/3 data point to 9 ohm-cm. (It is probably inappro- 
priate to compare for the CZ cell with the straight-line fit for higher - 
resistivity FZ material. ) 

4. 1.4 Electron Energy Dependence of Damage Coefficient 

Average damage coefficients for three resistivities are shown in Figure 14 

for the three electron energies utilized in this study. (Data labeled FZ 0. 1 

represent an average for both FZLD 0. 1 and FZHD 0. 1 samples, ) For 

comparison, findings of Downing et al. ^9 f or 3 ^ 3 ohm-cm solar cells are 

also presented. Upon comparing their data with the present FZ 2. 5 findings, 

it is seen that the solar -cell- to -bulk-material ratio tends to decrease 

with increasing electron energy. We also note that K-^ falls off more sharply 

with decreasing energy for the present data than for that of Downing, and thus 

our findings deviate more strongly than his from theoretical predictions^ of 

the energy dependence of K . 

L : 

4. 2 Results of 1 0 -Me V Proton Irradiations 

Data for 10-MeV proton-irradiated bulk samples and solar cells are listed 
in Table XI and shown in Figures 15 and 16. In Figure 15, saturation effects 
are again evident, particularly for low- resistivity samples (discussed in 
Section 4. 3). Additionally, a fluence dependence of the degradation rate is 
evident for solar cells. Also shown are data obtained at 9 x 10^ p/cm^ for 
two unilluminated bulk specimens. These data indicate slightly lower 
values for unilluminated samples, but more data are required before a 
more definitive statement can be made. 

11 7 

The first 10-MeV proton irradiation was performed at a fluence of ~ 10 p/cm 6 , 
followed by ~ 3 x 10 11 (cumulative) and 1.2 x lO 1 ^ p/cm . Carrier removal 
effects were negligible at this fluence for all resistivities examined. At this 
point, saturation effects were quite evident and thus it did not appear to 
be informative to irradiate the original group of bulk samples to an even 
higher fluence. However, it was obvious that data for fluences lower than 
10 * 1 p/cm^ would be beneficial for the purpose of obtaining damage co- 
efficients, particularly for low resistivity specimens. A subsequent 
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Figure 13. Damage coefficient vs. resistivity for 2. 5-MeV electron- irradiated 
specimens. 
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Diffusion length as a function of fluence for 10-MeV proton -irradiated bulk silicon specimens. 
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Figure 15. The quantity (L - Lq ) vs. fluence for 10-MeV proton- irradiated 
bulk silicon specimens and silicon solar cells. 
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Figure 16. Damage coefficient vs. resistivity for 10-MeV proton- irradiated 
bulk silicon specimens and silicon solar cells. 


56 


RESISTIVITY (ohm cm) 



irradiation was performed at ^3 x 10 p/cm . As seen in Table XI, the 
number of bulk samples irradiated at this fluence was small. Additional 
irradiations at relatively low proton fluences on a larger group of low- 
resistivity bulk samples would clearly have been useful in subsequent 
determinations. 

In Figure 16, is shown vs. resistivity. Values shown for low-resistivity 
bulk samples are for the two lowest fluences, and the effects of the onset of 
saturation are evident (lower K ^ for the higher fluence). A line with a slope 
of -0.44 provides a reasonable fit to several groups of data points/ Shown 
for comparison are data obtained by other workers. 2 /31-3 ^ i _MeV 

data of Anspaugh and Carter ^ * are somewhat lower than all other values. 
(The bars shown correspond to the range of values given in Reference 31. ) 
However, their measurements of diffusion length were made at a high in- 
jection level, and thus one expects a larger L, and hence a smaller K^, 
due to the injection-level dependence of diffusion length. Data of Meulenberg 
and Treble^® are also shown, and the bar shown corresponds to the fluence 

Q O 

dependence of Kj^ that they observed. Data obtained by Rosenzweig-^ agree 
well with the straight-line fit through several of the present solar-cell data 
points and also data of Meulenberg and Treble. Findings of Denny and 
Downing^ are also presented in Figure 16. 

Because of observed saturation effects for bulk samples and fluence-depen- 
dent K-^-values for solar cells, one must exercise caution when comparing 
device and bulk damage coefficients in Figure 16. At the lowest fluence, 
the solar-cell-to-bulk-material ratio for CZ 0. 1 is ~3, and this ratio 
for FZ 0. 1 is ~1. 7. Upon comparing solar-cell Kj^-values obtained at the 
lowest fluence with bulk material damage coefficients obtained from unity- 
slope fits to the data (Figure 15), the following ratios result: for FZ 2. 5, 

~ 4. 1; for CZ 2. 5, ~ 3.4; for FZ 13, ~ 7. 5. The latter value is anomalously 
high. As seen in Figure 16, K ^ for FZ 13/SC35 does not agree with other 
solar cell data presented; the reason for this discrepancy is unknown at 
present. 

Because of the paucity of low- resistivity data for bulk samples at the lowest 
fluence and because of saturation effects observed at higher fluences, no 
conclusion can be reached regarding the effect of dislocation density on K-^ 
for proton-irradiated, bulk specimens. However, for low- resistivity solar 
cells, data for CZLD 0. 1 units and CZHD 0. 1 devices are in good agree- 
ment (Figure 15), which suggests that dislocation density is an unimportant 
factor for 10-MeV proton damage. (A similar statement can be made for 
FZLD 0. 1/SC4 and FZHD 0. 1/SC2.) Regarding differences between FZ 
and CZ material, Figures 15 and 16 reveal once again that FZ solar cells 
are more radiation tolerant than their CZ counterparts. For example, 
the CZ-to-FZ ratio for 0. 1 ohm- cm cells at 1.2 x l(/2 p/cm^ is ~ 1. 5. 
For bulk material, upon comparing the K ^ value for CZ 6 in Figure 16 
with a straight-line fit to data for FZ 2. 5 and FZ 13, a ratio of ~ 1. 7 is 
obtained. 
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4. 3 Discussion of Saturation Effects 

Saturation effects were observed for irradiated bulk silicon 
a departure from linearity was noted in plots of the quantity 
versus fluence. Examination of Tables VIII through XI and Figures 8, 10, 12, 
and 15 reveals that the onset of saturation begins in all cases when L is de- 
graded to a value on the order of 30 pm. Continued irradiation reduces the 
apparent diffusion length somewhat to a lower limit of ~ 20 pm, and no further 
degradation is noted. Although such saturation behavior has been previously 
observed, ^ this effect is clearly not a manifestation of a saturation in 
the production of recombination centers because solar cells continued to 
degrade at high fluences while bulk samples saturated. 

We have had experimental indications that saturation is attributable 
to the effects of radiation-induced trapping centers. For heavily degraded 
samples, it was difficult to achieve a traps -filled condition using background 
lamps. That is, a slight decrease in observed SSPC signal was noted as 
the background lamp intensity was increased from about one-half of the 
maximum value up to the maximum. This was not the case for samples in 
the linear degradation region. Additionally, photoconductivity-decay life- 
time measurements were made for two CZ 2. 5 samples: one that had ex- 
hibited saturation and another that had not. Distinct differences in behavior 
were noted for the two samples. The saturation sample exhibited strong 
trapping behavior consistent with the presence of several types of centers. 
These centers possessed a range of trapping time constants, including times 
comparable to and also considerably longer than the lifetime. For this 
sample, because of trapping effects we were unable to extract a meaningful 
lifetime value from the observed decay even when a background lamp was 
employed. On the other hand, for the nonsaturation sample only slight 
trapping effects were noted. An exponential decay was readily extracted from 
the observed signal, and the lifetime thus obtained was in good agreement 
with that determined using the SSPC apparatus. 

The above observations suggest that in a heavily irradiated bulk specimen, 
i. e. , one in which L. has degraded to ~ 30 pm, traps introduced by radiation 
begin to have an effect on the SSPC signal. If trapping time constants are 
comparable to the lifetime, and the number of traps present is appreciable, 
then the SSPC method will not yield an accurate lifetime value. It will be 
erroneously high due to trapping, and background illumination does not 
provide a practical solution to this problem. For the present experiments, in 
most cases a region of linear degradation of L was observed at lower fluences, 
which enabled the determination of damage coefficients. Additional study 
would be necessary to gain a more detailed understanding of observed 
saturation effects. However, such an understanding is not required to achieve 
the technical goals of the present investigation. 

4. 4 Comparison of Steady-State Photoconductivity and Steady-State 
Surface Photovoltage Measurements 

A comparison was made of diffusion lengths measured by two techniques: 


samples in which 
(L-2 - L ” 2 ) 
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steady- state photoconductivity and steady- state surface photovoltage. 

Four irradiated bulk samples were examined. Following SSPC measure- 
ments, sample surfaces were prepared for the SPV technique in the 
following manner. Surfaces were lapped (6 |j,m) then polished with Lustrox 
1000. After this, surfaces were etched with CP4A followed by HF-H^O 
(1:4). Photon intensity vs. reciprocal wavelength is presented in Figure 
17 for the four bulk specimens. Extrapolated least-squares fits to the 
data intercept the abscissa at points numerically equal to bulk diffusion 
lengths. Table XII compares SSPC and SPV measurements for these 
specimens. 

TABLE XII 

Predicted SSPC L-values 
based on an extrapolated 


Sample 

SPV 

SSPC 

unity- slope fit 

FZLD 0. 1/1 

97. 4 Jim 

76. 5 ,um 

(not applicable) 

CZ 2. 5/3 

76. 4 

71. 2 

f 1 

FZ 13/18 

14. 4 

24. 7 

20. 8 jim 

CZHD 0. 1/8 

6. 2 

23. 4 

12. 0 


Comparison of diffusion lengths measured by SSPC and SPV techniques. 


It is seen that agreement between the two methods is reasonably good for 
long diffusion lengths but is rather poor at short diffusion lengths. How- 
ever, the two SSPC short diffusion-length values were measured on samples 
that^were heavily irradiated and exhibited saturation in plots of the quantity 
(L” - Lq~ 2 ) vs. fluence. As discussed above, such behavior is most likely 

due to trapping effects. In an attempt to make a more appropriate SSPC- 
to-SPV comparison, linear-degradation behavior observed at low fluences 
was extrapolated to the fluence of interest and an expected value of L (i. e. , 
in the absence of saturation) then calculated. These predicted L-values 
are also listed in Table XII. Agreement between SSPC and SPV values 
improves using this procedure, but significant differences still occur, 
particularly for CZHD 0. 1/8. 

In Sections 4. 1 and 4. 2, differences between damage coefficients for solar 
cells and bulk samples were noted. A comparison is made here of SPV 
and SSPC findings for the purpose of attempting to account for these 
differences. Figure 18 compares SSPC and SPV data for FZLD 0. 1/1 in 
a plot of the quantity (L - ^ - Lg“^) vs fluence. Also shown is a similar 
comparison for CZ 2. 5/13. It is seen that SSPC and SPV data agree quite 
well for CZ 2. 5/13. The solar cell data point for CZ 9/SC4 in Figure 
18 indicates the rather large device-vs -bulk discrepancy that needs to 
be accounted for. For sample FZLD 0. 1 / 1 , the SPV data point yields 
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Figure 17. Photon intensity (arbitrary units) vs. reciprocal absorption coefficient for 
irradiated bulk silicon specimens. Data shown were obtained by the steady- 
state surface photovoltage method for determining diffusion length. Dashed 

lines are least- squares fits, and the intercepts yield diffusion lengths. 
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a value even smaller than that for SSPC, and such behavior is in a direc- 
tion opposite to that required to account for differences between solar cells 
and bulk material. 


Figure 19 presents a similar comparison for the two short-diffusion-length 
samples of Table XII. For FZ 13/18, the SPV data point lies between the 
extrapolated SSPC value and that expected for a solar cell. To make this 
comparison, it is necessary to extrapolate linearly the behavior for FZ 13/SC3 5 
to a fluence of 1.2 x 10^ p/cm^. (A fluence-dependent damage coefficient 
would invalidate this comparison, however. ) For sample CZHD 0. 1/8, it 
is seen that the SPV result will give a value considerably larger than 
that expected for CZ 0. 1 solar cells at the fluence in question. 

It is encouraging that SSPC and SPV L.- values are in reasonable agreement 
for the two longer-diffusion-length cases examined. However, the dis- 
agreement noted for the two shorter-L samples indicates that additional 
comparisons of the two techniques in question, along with comparisons 
among other measurement approaches, will be required before observed 
differences between solar-cell and bulk-material damage coefficients can 
be resolved. 


4. 5 Annealing Study for Low-Resistivity Samples . 


An annealing study was performed for selected bulk silicon samples that 
had been irradiated at different energies to different fluences. One group 
of samples was annealed at room temperature following irradiation, and 
a second group was annealed at 60°C. Determinations of diffusion length 
using SSPC measurements were made for each sample at various times 
during the annealing period. Data were analyzed in terms of an unannealed 
fraction, defined through the relation 


Unannealed Fraction 


(,/L po S ,>- ,1/L 0> 


(15) 


In this expression, Lq is pre-irradiation diffusion length, Lp 0st is post- 
irradiation diffusion length measured after ^1 h at room temperature, 
and is post-ir radiation diffusion length measured at room temperature 
after a time t at the annealing temperature. Thus, a value of unity for the 
unannealed fraction occurs by definition at 1 h after irradiation. 


Room-temperature annealing data for twelve samples are shown in Figures- 
20 and 21. Very little annealing occurred during 500 h at room temperature 
for most of the samples. Sample FZHD 0. 1/5, however, exhibited a sig- 
nificant amount of recovery. This sample was irradiated with 0. 5-MeV 
electrons to a fluence of 1 x 10^ e/cm . It is interesting to note that 
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The quantity (L, - Lq ) vs. fluence for 10-MeV proton- and 1.0-MeV 

electron-irradiated bulk silicon specimens and silicon solar cells. A 
comparison is made of SSPC, SPY, and Co-60 data. 


Figure 19. 
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Figure 20. Unannealed fraction vs. time at room temperature for irradiated bulk 
silicon specimens. ' • ' 


64 


TIME (hours) 



CN 



oo o 

o o o 

NOI10VH3 Q3"IV3NNVNfl 


CN 

o 


o 


Figure 21. Unannealed fraction vs. time at room temperature for irradiated bulk 
silicon specimens. 
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the three other FZHD 0. 1 specimens studied exhibited more recovery at 
500 h than the other material types. Very little reverse annealing was 
observed. 

Figure 22 presents data for samples annealed at 60°C. After 500 h, six 
out of the seven samples have recovered somewhat, with an FZHD 0. 1 
sample once again exhibiting the most annealing. Reverse annealing was 
evident for one specimen (CZLD 0. 1/10). 

It should be noted that data in Figures 20 through 22 are for samples with 
differing irradiation histories. For certain samples, a plot of the quantity 
(L~ - L 0 ' 2 ) versus fluence would reveal that irradiations for these speci- 

mens were discontinued at a fluence in a region of linear degradation. For 
other samples, irradiation was discontinued in a saturation region and 
measured diffusion length was most likely influenced by trapping (see Section 
4. 3). In spite of these sample differences, with only one exception strong 
variations in annealing behavior among the eighteen samples of Figures 
2 0 through 22 are not apparent. 

4. 6 Solar Simulator Results for Irradiated Cells 

Figure 23 presents data showing the degradation with fluence of short- 

circuit current (I sc ), open-circuit voltage (V oc ), and maximum power 

output (P max ) for a 0. 5-MeV electron-irradiated 9 ohm-cm CZ solar 

cell. Data were obtained using a solar simulator, as discussed in Section 

2.1, and are normalized to pre -irradiation values (I gco , V oco , and P maxo ) 

of the cell parameters obtained from current- voltage characteristics. Pre- * 

irradiation cell properties are given in Table VII. The cell thickness for CZ 9/SC2 

was 0. 0285 cm. Downing's data 2 ^ indicate that a 1-MeV electron is 

"°3. 7 times as effective in degrading an n-on-p solar cell-as a 0. 5-MeV 

electron. Upon referring to the Solar Cell Radiation Handbook, ^ it is 

seen that for a 7-13 ohm-cm cell of the thickness employed here, I gc 

degrades to ^86% of its pre-irradiation value at a fluence of 8. 1 x 10 ^ 

1-MeV electrons /cm 2 (refer to Figure 3. 11 of Reference 11). At a fluence 
of 3 x 10^ 0, 5-MeV electrons /cm , which is 3. 7 times the 1-MeV fluence 
just given. Figure 2 3 reveals that I gc has degraded to 85% of its initial 
value. This agreement between data given in Reference 11 and the present 
findings suggests that the fluence measurements for 0. 5-MeV electron ir- 
radiations were accurate in the current investigation. (Good agreement 

is also found when the above calculation is performed for V and P . ) 

oc max 

Figure 24 presents data exhibiting the degradation of short-circuit current 
with fluence for 1-MeV electron-irradiated solar cells of several resis- 
tivities. Corresponding plots for open-circuit voltage and maximum power 
output are shown in Figures 25 and 26, respectively. The thickness of 
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Figure 22. Unannealed fraction vs. time at 60 C for irradiated bulk silicon specimens 
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Figure 23. Relative solar cell parameters vs. fluence for a 0. 5-MeV electron- 
irradiated unit. 
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Figure 25. The quantity V /V vs. fluence for 1. 0-MeV electron- irradiated silicon 
, oc oco 

solar cells. 
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Figure 26. 


The quantity P max / p maxo vs. fluence for 1.0-MeV electron-irradiated 
silicon solar cells. 
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device CZ 9/SC3 is 0. 0307 cm, and by comparing the present findings with 
expected behavior for cells irradiated with 1-MeV electrons, ^ good agree- 
ment is found, which once again suggests that fluence measurements were 
accurate. Figure 24 reveals that I degraded to ' v, 62% of its pre-irradiation 
value at 3 x 10^ e/cm^ for 0. 1 ohm-cm cells, as compared to ^75% for 
higher-resistivity units. For low-resistivity devices, no dependence of 
I degradation on dislocation density is apparent. 


Figure 25 shows that V q degraded to ~ 88% of its initial value for CZ 2. 5 and 
CZ 9 cells at 3 x 10 15 e [cm^, whereas the corresponding value for low-resis- 
tivity cells is — 91%. In Figure 26, it is seen that in terms of normalized 
maximum power output higher resistivity cells are superior. For such units, 

at 3 x 10l 5 e/cm^ / P = 0. 65, whereas the low resistivity value is 

max maxo 

^ 0. 53. (The one CZLD 0. 1 cell tested appears to be slightly superior to the 
other 0. 1 ohm-cm cells; however, in terms of absolute power output, this unit 
is actually inferior to the others.) Thus, solar simulator findings are con-.' 
sistent with damage coefficient determinations reported above in that low- . * 
resistivity solar cells are inferior to higher resistivity units in terms of 
radiation tolerance. 


Figure 27 presents degradation of solar cell parameters with fluence for 
a 2. 5~MeV electron-irradiated 9 ohm-cm device with a thickness of 0. 0302 
cm. Downing's data^ ^ indicate that a 2. 5-MeV electron is ^4. 1 times as 
effective in degrading an n-on-p cell as a 1. 0-MeV electron. Upon taking 
this difference into account and then comparing the data of Figure 27 with 
that in the Radiation Handbook^ in a manner similar to that employed 
above for 0. 5-MeV electrons, good agreement is found. This agreement 
again suggests accurate dosimetry for the present electron irradiations. 


Figure 28 presents the degradation of short-circuit current with fluence 
for 10-MeV proton-irradiated solar cells of two resistivities (0. 1 and 
2. 5 ohm-cm). Corresponding plots for open-circuit voltage and maximum 
power output are shown in Figures 29 and 30, respectively. As a check 
on fluence accuracy, consider data for 2. 5 ohm-cm cells. Data in the 
Radiation Handbook^ (see Figure 3.6 of this reference) indicate that a 
10-MeV proton fluence of 3 x 10 ^ p/crn^ will degrade I gc by the same 
amount as a 1-MeV electron fluence of 2 x 10^ e/cm^. From Figure 
28, we find that I sc for CZ 2. 5 cells degrades to ^ 0 . 79 of its initial 
value at 3 x 10^ p/cm^. From Figure 24, interpolation yields ~0. 76 
for this same quantity at 2 x 10^ e/cm^ for a 1. 0-MeV electron-irradi- 
ated CZ 2. 5 cell. This agreement suggests that fluence determinations 
for 10-MeV proton irradiations were accurate. 


Figure 28 indicates that normalized I degrades more severely for 0. 1 ohm-cm 
cells than for higher - resistivity units. However, as shown in Figure 29, 
the opposite is true for V degradation. In Figure 30, it is seen that 
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Figure 27. 
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Figure 28. The quantity I /I vs. fluence for 10-MeV proton- irradiated silicon 
solar cells. SC SC ° 
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Figure 30. The quantity P /P vs. fluence for 10-MeV proton-irradiated 
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silicon solar cells. 
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for three out of the four low-resistivity cells examined, P /p was 

_ , , max' maxo 

somewhat less than for 2.5 ohm-cm cells which is consistent with expecta- 
tions based on determinations. For CZLD 0. 1/SC3 this quantity was 
identical to that for CZ 2.5/SC3. However, in terms of absolute quantities 
(use Table VII in conjunction with Figure 30), the power output for CZLD 
0. 1/SC3 was the lowest of all six cells shown in Figure 30. (Caution must be 
exercised in interpreting the normalized data of Figures 23 through 30. ) 

After 1.2 x 10^ protons/cm^, none of the four low resistivity cells had a 
maximum power output as large as that for the two CZ 2.5 devices. 
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SECTION 5. 0 


DISCUSSION AND CONCLUSIONS 

One finding of practical importance in the present investigation is the observa- 
tion of damage coefficients that increase with decreasing resistivity (p) over a 
broad range of resistivities. The employment of low-resistivity solar cells in 
a space radiation environment therefore may not be advisable. However, whether 
low-resistivity cells may or may not be advantageously utilized in space applica- 
tions will depend on both the amount of radiation-induced degradation anticipated 
and the pre- irradiation properties of such devices compared to those for higher- 
resistivity units. The obtaining of a basic mechanistic understanding of the 
resistivity dependence of is an intriguing question that has been addressed 
by others. 29,36 -^e ^ave a i so attempted to understand such behavior, and our 
preliminary analysis is now given. 

For all three electron energies examined, empirical fits to the data yield the 
following approximate relation: 

K l « p" 2/3 . (16) 

Hole density p Q was determined for p-type material over the resistivity range 
from 0. 1 to 10 ohm-cm. Additionally, damage coefficient is proportional to 
(L - ^ - Lq”^), which can be interpreted as the reciprocal of radiation -induced 
carrier lifetime t . Thus, Kj^ « 1/ T. Equation 16 was then employed to 
obtain a plot of t vs P q > which is shown as the solid curve in Figure 31. This 
curve is thus representative of the experimental data obtained here for three 
electron energies. 

We next attempted to fit the solid curve using the Hall-Shockley-Read (HSR) 
model. The HSR lifetime expression for a single recombination level in p- 
type material at low injection levels is 


T 


p o + P 1 

C „ Np o 


n + n, 
o 1 


c Np 
p o 


(17) 


where all terms have their usual meaning. As p Q increases, lifetime decreases 
due to the increased recombination rate. For ease of analysis and for purposes 
of illustration, a recombination level in the lower half of the bandgap was as- 
sumed. For any reasonable capture -probality- ratio assumed. Equation 17 then 
simplifies to 


T 



(18) 


We attempted to fit Equation 18 to the solid curve in Figure 31 by treating the 
multiplicative term l/c n N as an adjustable scale factor. A reasonable fit was 
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RADIATION-INDUCED LIFETIME (relative units) 



Figure 31. Radiation-induced lifetime vs. hole density. It is demonstrated that 

recombination through two Hall-Shockley-Read levels can be employed 
to qualitatively explain the observed dependence of damage coefficient 
on resistivity for electron -irradiated boron-doped silicon. The re- 
combination levels given here are not unique to the experimental data 
and are presented for purposes of illustration only. 




obtained for a level at E^ - E v — 0. 12 eV. However, an exact fit can be ob- 
tained by using a two-level HSR model, as shown in the figure. By assuming 
that each level can be expressed in the form of Equation 18 and by adding the 
contributions of each level reciprocally to obtain t, the solid curve represent- 
ing the experimental data can be fit exactly. As shown in Figure 31, the level 
more effective for higher resistivities is at E r - E v a: 0. 18 eV, whereas the 
level more effective at lower resistivities is at E r - E v s: 0. 08 eV. The ratio 
of the multiplicative scale factors for these two levels is ~ 22, with c^N for the 

0. 18 eV level being larger. 

No physical significance should be attached to the above recombination param- 
eters. Levels in the upper half of the bandgap could be effective and/or three 
or more levels may be effective rather than only two. (It should be noted that 
for the two-level model employed here, a satisfactory fit to the data requires 
a shallow level close to either the valence band or the conduction band. ) Our 
purpose here was only to illustrate that a two-level HSR model describes the 
present data quite well. In order to obtain meaningful parameters from such a 
model, the temperature and injection-level dependences of radiation- induced 
lifetime are required. 38 Wilsey^ has examined the temperature dependence 
of K l for electron-irradiated p-type silicon, and it is clear that a one-level 
model cannot account for his data. This statement can also be made for findings 
of Faith. in terms of gaining increased understanding of the dependence of 
on resistivity, it would be useful to measure the temperature and injection- 
level dependences of damage coefficient for samples of several resistivities 
and then attempt to fit the data with two (or more) HSR levels. 

Several of the major observations made and conclusions reached in this in- 
vestigation are the following: 

1. Diffusion -length damage coefficients increase with decreasing 
resistivity for boron-doped silicon. For low-resistivity solar 
cells, this decrease in radiation tolerance must be weighed 
against potential advantages when considering such devices for 
utilization in a space radiation environment. 

2. For 0.5-, 1.0-, and 2. 5-MeV electron bombardment, empirical 
fits to experimental data can be approximately expressed as 

K l <x p _ 2/3 f or o. 1 ^ 20 ohm- cm. F or 10-MeV proton 

bombardment, an empirical fit of the form « p“ ^ was 
found to describe the data reasonably well. 

3. The dependence of damage coefficient on resistivity can be 
qualitatively accounted for quite well using a two-level Hall- 
Shockley-Read model. 

4. Damage coefficients for solar cells were observed to be larger 
than their bulk-material counterparts. 
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Bulk samples and solar cells prepared from float-zone material 
were generally observed to be more radiation tolerant than their 
Czochralski counterparts at all resistivities examined. 

6. No dependence of damage coefficient on dislocation density was 
apparent for 0. 1 ohm-cm bulk samples and solar cells. 


Regarding item 4, during the final stages of preparing this report a comparison 
was made of measured pre-irradiation diffusion lengths for bulk samples and 
solar cells fabricated from the same ingot. Bulk values were found to be a factor 
of 1. 7 ± 0. 3 larger than their device counterparts. This value is an average 
obtained by comparing 110 bulk samples and 33 solar cells fabricated from 
7 ingots. Referring to Equation (1), it is seen that if pre- and post-irradiation 
values of diffusion length for a given bulk sample were a factor of 1. 7 larger 
than the corresponding solar-cell values, then the solar cell value would 
be a factor of 2. 9 larger than the bulk damage coefficient (assuming equal 
fluences). This latter factor is comparable to solar-cell-to-bulk-material 
damage- coefficient ratios experimentally observed in this study, which suggests 
that a systematic error is responsible for the device-bulk differences noted. 

Several items that should be considered when attempting to quantitatively 
account for the observed discrepancy are: a) surface photovoltage diffusion- 
length measurements were found to be in reasonable agreement with SSPC 
measurements on two bulk samples with relatively long diffusion lengths (refer 
to Section 4.4); b) the present 1.0-MeV solar-cell damage coefficients are in 
reasonable agreement with those obtained by other workers (refer to Figure 11); 
c) the fluence dependence of for solar cells noted here and also by other 
workers is in a direction to bring device and bulk damage coefficients into 
agreement at high fluences; (this statement assumes that bulk values are 
not fluence dependent); d) accurate determination of the generation rate in a 
penetrating -radiation diffusion -length measurement for solar cells is not 
straightforward. We presently suspect that previous workers and ourselves 
have not properly accounted for all of the important factors in penetrating 
radiation solar-cell measurements. Some of the difficulties associated with 
determining the generation rate for a gamma- irradiation diffusion -length measure- 
ment have been previously considered. ' ^ It seems clear that additional 
work will be required to obtain a satisfactory explanation for observed device- 
bulk differences. 
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APPENDIX 


STATISTICAL ANALYSIS OF DATA" 

Electron Damage Coefficients 

The present study has examined radiation effects on four types of 0. 1 ohm- 
cm boron-doped silicon ingots, including the four possible combinations of 
high and low oxygen content and high and low dislocation density. For 
electron-irradiated specimens, Figures 9, 11, and 13 indicate that there 
is no significant dependence of damage coefficient on dislocation density. 
However, a dependence of Kj^ on growth method (i. e . , oxygen content) is 
evident. In this section, an attempt is made to quantitatively examine the 
statistical significance of the dependence of electron damage coefficient on 
oxygen content (OC), dislocation density (DD), and electron energy (EE). 
Additionally, the dependence of pre -irradiation diffusion length on oxygen 
content and dislocation density is examined. 

Statistical analyses were performed by utilizing analysis of variance pro- 
cedures outlined by Hicks. ^0 por each 0. 1 ohm-cm ingot type at a given 
electron energy, values of for three individual samples were employed. 
Table XIII lists damage coefficients used in the present analysis. This 
table is arranged in a matrix corresponding to the four possible material 
combinations. (We associate a low oxygen concentration with float- zone 
material and a high oxygen concentration with Czochralski material. ) Three 
variables, or factors, are of interest here: oxygen content, dislocation 
density, and electron energy. The first two factors were examined experi- 
mentally at two levels and the last at three levels. This situation is termed 
a 2 x 2 x 3 factorial experiment with three replications per cell. 

The damage coefficients in Table XIH are summarized in a format usually 
associated with factorial designs. Because of the wide range of values 
observed (almost two orders of magnitude) and observation of the fact that 
the standard deviation of the three numbers in each box is proportional to 
the absolute value of the numbers, the logarithms of the damage coefficients 
were calculated. These are also listed in Table XIII. A second reason for 
this approach is that the analysis of variance statistical procedure requires 
that the error variance of each observation be equal. This is much more 
true of the logarithms than of the actual values. 


Analysis presented in this Section was performed by S. M. Sidik and 
C. R. Baraona of NASA-Lewis Research Center. 
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FZLD 


FZHD 


CZLD 


CZHD 








4. 44 

4. 75 

14. 0 

11. 9 


(0.65) 

(0. 68) 

(1. 15) 

rasa 

0. 5 MeV 

4. 13 
(0. 62) 

4. 20 
(0. 62) 

14. 0 
(1. 15) 

11.2 
(1. 05) 


4. 00 

4. 20 

13. 3 

10.2 


(0.60) 

(0. 62) 

(1. 12) 

B&SII 


42. 8 

45.6 

72. 0 

85. 0 


(1.63) 

(1. 66) 

(1.86) 

(1. 93) 

1. 0 MeV 

42. 8 

38.2 

72. 0 

85. 0 

(1.63) 

(1. 58) 

(1. 86) 

(1. 93) 


42. 8 

38. 2 

67. 0 

67. 0 


(1.63) 

(1. 58) 

(1. 83) 

(1. 83) 


194 

196 

220 

237 


(2. 29) 

(2.29) 

(2. 34) 

(2. 37) 

2. 5 MeV 

173 
(2. 24) 

188 

(2.27) 

261 
(2. 42) 

228 

(2.36) 


176 

200 

241 

249 


(2.46) 

(2. 30) 

(2. 38) 

(2. 40) 


TABLE XIH 

Damage coefficients obtained at three electron energies for 
samples from four 0. 1 ohm-cm B-doped Si ingots. All damage 
coefficients listed are in units of electrons - ^ and should be 
multiplied by lO - ^. Numbers shown in parentheses are the 
logarithms of the listed damage coefficient values. 
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Table XIV presents results of the analysis of variance for the data (logarithms) 
of Table XIII. Listed are the sources of variability in the experiment, the 
degrees of freedom, the observed mean squares, the expected mean squares, 
the calculated F ratio test statistic, the standard table value of the F statistic 
(95% confidence level), and the resulting significance of each source of variance 
in the experiment. 

Table XIV is divided into three parts by horizontal lines. This emphasizes 
the structure of the expected mean squares, because these structures 
determine which F tests should be used. For instance, the expected mean 
squares for EE, OC x EE, DD x EE, and OC x DD x EE are all of the form 
dg plus one other source of variability, where is the error variance. 

Their significance may thus be tested by dividing those mean squares by the 
error mean square. These F statistics indicate that EE and OC x EE are 
highly significant (values far outside of most tables) and that DD x EE and 
OC x DD x EE are marginally significant at the 95% confidence level. 

The part of the table corresponding to OC, DD, and OC x DD is somewhat 
different. The expected mean squares are all of the form dg + 9 d| plus 
one other source of variability. The symbol dg denotes the variance of 
damage coefficients resulting merely from using different ingots. These 
mean squares should be tested by dividing by a mean square which has an 
expected value of dg + 9dg , but there is no such estimator available in the 
present experiment. The only alternative available is to assume there is 
no interaction between DD and OC. In this event, OC and DD may be tested 
for significance by dividing the mean squares by the mean square for the 
DD x OC interaction. This was done and shows that OC is highly significant 
and that DD is highly insignificant. In Figure 32, average values of damage 
coefficients in Table XIII are plotted versus electron energy, and the trend 
of the data graphically verifies the statistical results. 


Pre -Irradiation Diffusion Length 

For the same four ingots discussed above, pre-irradiation diffusion lengths 
for 18 samples of each were statistically examined in a 2 x 2 factorial 
experiment. (See Table VI for numerical values of diffusion length.') As 
in the damage coefficient analysis, because of the manner in which the 
data were obtained there is no direct test for the effect of OC, DD, or 
their interaction. If the interaction is assumed not to exist, then DD and 
OC may be tested for significance by computing the ratio of their sums 
of squares and the interaction sum of squares. This was done and it was 
found that neither DD nor OC have a significant effect. This result may 
be due either to the fact that OC and DD do not, in fact, have any significant 
effect, or to the possibility that there is indeed an interaction between the 
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listed in Table XIII. 



DAMAGE COE FFECIENT (electrons 


A — - FZLD 0.1 
O — - FZHD 0.1 

O CZLD 0.1 

□ CZHD 0.1 

DATA POINTS ARE 
AVERAGE VALUES 


0 1.5 

ELECTRON ENERGY (MeV) 

Figure 32. Average damage coefficient vs electron energy for four 

0. 1 ohm-cm B-doped Si ingots. Average values are based 
on the data listed in Table XIII. 





two factors. If the latter were the case, then these tests would not be as 
sensitive as they should be. Therefore, no definite conclusion can be made 
concerning these factors. 

One conclusion can be drawn, however. The ratios of the mean squares 
due to OC, DD, and OC x DD were calculated and all were found to be 
highly significant. This indicates that there are significant differences 
among the four ingots with regard to pre -irradiation diffusion length (which 
is obvious from the data in Table VI). This occurs because the error mean 
square estimates the variability of diffusion length among samples within 
an ingot while the mean squares due to OC, DD, and OC x DD estimate the 
variability among different ingots. It is not possible to state the cause of 
the difference on the basis of the present experiments. 

Pre -irradiation data in Table VI suggest, on the surface, that oxygen 
content has a strong effect on diffusion length. However, past experience 
indicates that this is not necessarily the case. Lifetime (or diffusion 
length) is determined by the number and type of recombination centers 
present. In order to test properly for the dependence of pre -irradiation 
lifetime on oxygen content and dislocation density, all other properties 
for the ingots examined would have to be the same. Two factors that can 
affect lifetime include the impurities present in a given ingot (other than 
dopant concentration) and heat treatment history. For the present four- 
ingot comparison, one would have to characterize each in considerable 
detail and determine whatproperti.es are different other than dislocation 
density and oxygen content. Additionally, in an ideal experiment one would 
have to study a large number of ingots. In the absence of a detailed ingot 
characterization and because of the limited number of ingots employed, 
one cannot make a definite conclusion regarding the effect of oxygen or 
dislocation density on pre -irradiation diffusion length based on the present 
study. 
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